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OPERATOR SEPARATION OF VARIABLES FOR 
ADIABATIC PROBLEMS IN QUANTUM AND 

WAVE MECHANICS 

Belov V.V.* Dobrokhotov S.Yu.^ Tudorovskiy T.Ya3 


Abstract 

We study linear problems of mathematical physics in which the adiabatic approxi¬ 
mation is used in the wide sense. Using the idea that all these problems can be treated 
as problems with operator-valued symbol, we propose a general regular scheme of adia¬ 
batic approximation based on operator methods. This scheme is a generalization of the 
Born-Oppenheimer and Maslov methods, the Peierls snbstitution, etc. The approach 
proposed in this paper allows one to obtain “effective” reduced equations for a wide 
class of states inside terms (i.e., inside modes, subregions of dimensional quantization, 
etc.) with the possible degeneration taken into account. Next, applying the asymp¬ 
totic methods in particular, the semiclassical approximation method, to the reduced 
equation, one can classify the states corresponding to a distinguished term (effective 
Hamiltonian). We show that the adiabatic effective Hamiltonian and the semiclassical 
Hamiltonian can be different, which results in the appearance of “nonstandard charac¬ 
teristics” while one passes to classical mechanics. This approach is used to construct 
solutions of several problems in wave and quantum mechanics, in particular, problems 
in molecular physics, solid-state physics, nanophysics, hydrodynamics. 

Keywords: adiabatic and semiclassical approximation, Born-Oppenheimer method, Maslov 
operator methods. 

1 Introduction 

Many linear problems of mathematical and theoretical physics contain different spatio- 
temporal scales. Among them there are problems of molecular physics, problems concerning 
electron waves in crystals, wave propagation in media with rapidly varying characteristics, 
surface and internal waves in fluids, electron-phonon interaction, electromagnetic waves and 
quantum particles propagation in waveguides, etc. The main instruments for investigating 
such problems are contained in the adiabatic approximation, which is based on the idea of 
separation of “fast” and “slow” modes by means of “freezing” the slowly varying variables. 
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For instance, the “slow” variables describe the nuclear motion in molecules and the “fast” 
variables concern the electron motion, or the “slow” variables describe the longitudinal mo¬ 
tion and the “fast” variables describe the transverse motion in thin waveguides. Needless to 
say that there are many different versions of adiabatic approximation and thousands papers 
and monographs related to this approach and its applications in different fields of mechan¬ 
ics and physics. Among them we mention niiaisiiiiiniiniiziiHiEiiini Nevertheless, we 
take the liberty to present a general regular scheme of the adiabatic approximation sug¬ 
gested in mmaiini and to combine different approaches including the Born-Oppenheimer 
method, the Maslov operator method, the Peierls substitution, etc. From the mathemat¬ 
ical viewpoint, the equations (or the system of equations) describing all these phenomena 
have the same structure. Namely, following [Hj, these equations (systems) can be treated 
as equations with “operator-valued symbol.” Our idea of the study of such equations is 
not new: the asymptotic analysis of the original problem can be divided into two parts: 
(1) the “operator” reduction to simpler differential or pseudodifferential equations with the 
principal symbol, known in different fields of physics as an effective Hamiltonian, or therm, 
or dispersion relation, or mode, etc., and with corrections to this symbol; (2) asymptotic 
constructions of the solutions of this simpler reduced equation based on different variants 
of the semiclassical approaches, like the WKB-method, Born-method, oscillatory approxi¬ 
mations, ray expansions, the Maslov canonical operator, averaging, etc. Here we present 
the hrst part of this concept in the form of a regular rigorous algorithm (in §3), based on 
operator methods na. 

The result of the hrst step is the reduced equation; it has different names in different helds 
of physics, we call it the effective equation of adiabatic motion. We illustrate the “operator” 
reduction or the “operator separation” of variables by using the above-mentioned problems 
from different helds of physics and, in §4, present the corresponding equations for the wave 
functions of adiabatic motion. The examples given in items 4.1-4.4 of the paragraph were 
studied long ago, whereas the results of items 4.5-4.6 (as well as §5) were obtained by authors 
recently. 

To realize the second step, it is necessary to take into account that usually the original 
problem includes several parameters: some of them, like the transverse and longitudinal char¬ 
acteristic sizes of a waveguide or the ratio between the masses of light and heavy particles, 
allow one to use the adiabatic approximation and do not crucially correlate with the energy 
of adiabatic motion, and the other ones, like magnitudes of the external electromagnetic 
held, the momentum of the incoming wave in the scattering problem, etc., determine the 
energy. This fact implies diherent forms of (asymptotic) wave functions of adiabatic motion 
and, as a consequence, a redehnition of the principal symbol and the ehective Hamiltoni¬ 
ans depending on the relations between the above-mentioned parameters. In turn, it gives 
diherent types of characteristics (trajectories of Hamiltonian systems) which V. Maslov 
called “nonstandard characteristics” and which must be used in asymptotic constructions. 
We discuss the possible classihcation of these characteristics using, as the main example, 
the quantum wave propagation in thin (or nano) tubes. In spite of the fact that the given 
arguments seem to be natural, and in some way contain in physical literature, we didn’t hnd 
their systematic consideration. The methods for constructing asymptotic or exact solutions 
of the “redehned” equation for the wave functions of adiabatic motion are well known and 
here the results must be connected with a concrete physical problem. Therefore, we do not 
construct asymptotic solutions for most of the derived reduced equations and in §5 only 
briefly describe different solutions for the equations of quantum particles in nanostructures. 

The main results of the “operator separation of variables” are realized in formulas (ESI), 
(EH) and Eq. (ITTTD . Although they simply develop the approach of Born-Oppenhiemer, 
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Peierls and Maslov, nevertheless, they allow ns to consider a wide range of adiabatic prob¬ 
lems uniformly and in a rather compact form. We believe that this approach is very useful in 
different situations, since it gives not only a general regular scheme for deriving the reduced 
equations exactly but allows one to obtain qualitative and quantitative estimates of the range 
of applicability of any approximation. Naturally, the argument resulting in formulas (ESI) 
and (EZZD, the classihcation of different approximations, the relations between the adiabatic 
and semiclassical asymptotics, etc. can be better illustrated with a simple example. These 
considerations, some of which are well known in physics and some of which are well known 
in mathematics, are given in §2 and §§3.1-3.2. In §2, we present a minimal amount of the re¬ 
quired information from the operator calculus of noncommuting operators m (see also Hg). 
We point out that the facts from mi used here are not simply arguments of “mathematical 
justihcation and verihcation” type, but are completely constructive and developed algorithms 
well adjusted to the problems studied here. 

Finally we can formulate the main result of the paper in the following way. We suggest 
the regular asymptotic (adiabatic) procedure which allows one 1) to determine correctly the 
leading part of asymptotic solution corresponding to wide diapason of energies (or frequen¬ 
cies) 2) to construct and estimate if necessary the “adiabatic” corrections. Needless to say 
that in this work we understand an asymptotic solution in the formal sense, i.e. in the 
sense of small “right hand side” (discrepancy). The proof of the fact that the constructed 
asymptotic solution approximate some exact solution of original equation is out of scope of 
our consideration here and we touch this problem only very briefly. 


2 Differential and pseudo differential operators with a 
parameter and their symbols. Elementary formulas 
from calculus of noncommuting operators. 

We want to study some asymptotic solutions of (systems of) partial differential equations 
with small parameter ^ in the conhguration space with coordinates x = (a:i,... ,xp^) which 
can be written in general form: 


iTt = TiT. 


( 2 . 1 ) 


Here T {x, t) can be a scalar or vector function, is a partial differential scalar or ma¬ 
trix operator. It is convenient for us to present the operator 7Y as a function of noncom¬ 
muting operators —id/dx = {—id/dxi,...,—id/dxN) and x = {xi,... ,xn) and, gener¬ 
ally speaking, of time t: H = H{—id/dxi ,..., —id/dxN,xi,... ,XN,t), where the function 
7-f(pi,... ,pn, xi,..., xn, t) is usually called the symbol of the operator H. Actually, we shall 
consider the situation in which the function Ti can depend on the parameter p and also on 
some other ones. Quite often a small parameter p appears as a factor before the derivatives 
d/dxj, say, before d/dx \,..., djdxn-, n < N. It follows from the considerations given below 
that there is always a parameter p before d/dt. We denote the other variables hj yi,, ym, 
m = N — n. So hnally Eq. (EH) takes the form 


i/iTi = H I —iy 


A 

“dxi ’ 




d 


d 


, Xi, . . . , Xm i „ , 

OXn uyi 


—I 


A. 

dyn 


'i yij ■ ■ ■! ymt h) ^ (2-2) 


As the operators d/dxj and Xj, as well as operators d/dyk and yk do not commute, one 
has to agree about the order of action of Xj and d/dxj and, analogously, of yk and d/dyk. 
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The theory of functions of noncommuting operators is very well developed na, see also 
Ha 12111221 • For the completeness of consideration, let us present a minimal amount of the 
required information from the operator calculus and recall the terminology. 

First, let R{x,p) = Yfk=o^k {x)p^ be a polynomial in variables p with coefficients smooth 
in X. This function generates the operator R = Y^k=QRk{x){—ipL-^Y. The function R{x,p) 
is called the symbol of the differential operator R with a parameter p. It is clear that the way 
of the construction of the operator R by means of the symbol R is not unique. For example, 
one can build an operator R' = ^k{x)^ different from R. Using the Feynman 

1 2 2 1 

notation, we can write R = R{p,x) and R = R'{p,x), where the numbers above p and x 
determine the order of their action. (About other ways of ordering, for example, by Weyl, see 
jl4ji. In this work, we will always use the first way of ordering. Under this agreement, the 
definition of the differential operator is equivalent to the definition of its symbol. By letting 
the order of the polynomial k tend to infinity, one can obtain, at least, “naive” operators 
whose symbols are not polynomials. Such operators are called pseudo differential. Their 
rigorous definition is given by means of the “p-Fourier transform” PEj: 


A{p,x)^{x) = F^^^[A{p,x)[F^^m{x)]{p)]{x)- 


where the direct and inverse 
equalities: 

[F^^p{x)]{p) = 


“/i-Fourier transforms” 


F^^p and F^_ 


(27ri/i)"'/2 


e [F^^^p{p)]{x) = 


-27iip)'^R 


(2.3) 

are defined by the 




From now on, (,) is the inner product in the Euclidean space of the corresponding dimension. 

The replacement of operators by their symbols turns out to be very useful in practical 
calculations. As a result, the calculations concerning operators are replaced by significantly 
simpler work (which can be algorithmized) with symbols, i.e., functions (“with c-numbers”). 
Since, in asymptotic approaches, defining an operator is practically equivalent to defining 
its symbol, in the process of obtaining asymptotic formulas one can manipulate only with 
symbols and “recall” the operators corresponding to these symbols only in studying refined 
problems such as, for example, justification of the asymptotic accuracy of the solutions con¬ 
structed. Of course, the main difficulties in dealing with functions of operators arise due to 
the fact that the operators p and x do not commute. On the other hand, their commutator 
is ip, and it is small, which allows one to use asymptotic expansions in the constructions. 
In view of this fact, it is natural to consider the symbols R depending on the parameter p 
and to assume that R{p,x,p) = Ro{p,x) + pRi{p,x) + ■ ■ ■. Moreover, the right-hand side 
in this relation is understood as an asymptotic expansion in the parameter p. The func¬ 
tion Ro{p,x) is called the leading symbol or, sometimes, a Hamiltonian, and Rj are called 
jth-order corrections. 

The next generalization consists in the assumption that the symbol of the operator R 
may be an operator. A simple example appears in the situation in which R{p, x, p) is a 
matrix (or an operator acting in a finite-dimensional space). 

Example 1. The Klein—Gordon equation. Consider, for instance, the Klein-Gordon 
equation — P^'pxx + v{x)(p = 0 written in vector form for the vector function T = 

tA ^ / (p Y 
^2) VWv 


ip^t = 'H{—ipd/dx,x)^, 
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ip'^u = ^ 2 , 

ill'll!2t ^ 


( 2 . 4 ) 




The symbol of the operator T-l{p,x) is the 2x2 matrix function 


n{p,x) 


( ° 

\p‘^ + v{x) 



From this viewpoint, one can consider many fundamental physical equations like the 
Dirac and Pauli equations, the Lame equation in linear elasticity theory, the linearized 
hydrodynamics equations, etc. (If they include a small parameter in an appropriate way.) 
The appearance of a small parameter p before the derivative d/dx is very important in our 
constructions. As we have mentioned, there exist many problems with different scales in 
which a small parameter appears only in front of some derivatives. Problems of such types 
give the majority of nontrivial /i-differential operators with operator-valued symbols. 

Example 2. Molecular physics. Consider, for instance, the Schrodinger equation for 
two groups of particles: heavy atomic nuclei with mass M and light electrons with mass m. 
We denote the coordinates of nuclei and electrons by x' and y', respectively. Let us assume 
that Iq is the linear size of a molecule and do is the amplitude of nuclear oscillations. Thus 
the characteristic magnitude of the electron energy is Se ~ h^/{2mlQ). By physical reasons 
(stated by Born and Oppenheimer), the motion of a nucleus could be considered in the 
oscillatory approximation and its energy is /{2Md‘l) ~ fc(ig/2 with the elasticity 

coefficient k. To estimate k, one has to remember that, in the adiabatic approximation, 
the potential energy of a nucleus is the total energy of electrons , so k ^ d'^Ee/dx'^ ~ 
h^/(m/g). Thus we have k?/{2Md‘^) ~ kid'll[2x01^). From this, we obtain do/lo ~ (m/M)^A_ 
Oscillatory energies of nuclei and electrons relate as SnlEe ~ {m/M){ll/dl) ~ Let 

us introduce the parameter p = {rn/My/'^ and divide both sides of the Schrodinger equation 
by k?/{2mll). After passage to dimensionless variables x = x'/Iq, y = y' /Iq, the stationary 
Schrodinger equation takes the form 


TfT = ET, n = 


+ vix, 2/) ) ^ = ET. 


The symbol Td of the /x-differentiaP operator Td is again the operator 


W(p,i) = Y - Yy + v(x,y). 


(2.5) 


( 2 . 6 ) 


Usually, X are called slow variables and y are called fast ones. Ideologically close approaches 
to the determination of electron states [terms) in molecule one can hnd in jHl IHl IHl IT7j. 

Example 3. Quantum 2-D waveguide. One can meet an equation with closed struc¬ 
ture considering a “narrow” straight quantum waveguide. The word “narrow” means that 
the characteristic width of the waveguide do is much smaller than its length Iq. We introduce 
the small parameter p = do/lo- The dynamics of a spinless quantum (quasi)particle in a 
plane waveguide is determined by the 2-D Shrodinger equation with the potential v = v[x, y) 
inside the waveguide. Due to two different scales, there appear two different characteristic 
energies: the characteristic energy of the lower transverse levels (which is usually called the 
characteristic energy of the “transverse quantization”) and the characteristic longitudinal 
energy ey. One can estimate e_L from the uncertainty principle, which gives e± = fif/[2md/^. 
Let us introduce dimensionless variables x' = x/lo-, y' = y/do, t = {puj±)~^, uj±_ = e_\_/h and 

^Born and Oppenheimer in their famous paper pQ used the parameter x = U/I, which is the ratio 
do/lo ~ X of the characteristic wavelength to the wave function and the linear size of the molecule. 
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dimensionless potentials v' = v/e±. Then the corresponding Schrodinger equation takes the 
form (we omit the primes of the dimensionless variables): 

The symbol of the /i-differential operator is the operator (I2.t)|l with Ay = d‘^/dy‘^. 

Generalizations of the plane quantum waveguide are quantum thin tubes (nanotubes) 
and thin dims (nanohlms), and their symbols are matrix-operators if one includes spin into 
consideration. These more complicated examples, as well as several examples from other 
fields, will be considered later. 

Let us again stress that the definition of the symbol of a y-differential operator differs 
from the standard definition of the symbol of an operator without parameter. Namely, we 
construct the symbol taking only the slow variables into account. This is why the symbols of 
the /i-differential operator Td in Examples 2 and 3 are again differential operators acting in 
some appropriate Hilbert space with coordinates y. The transition to /i-differential symbols 
is a formalization of the idea of “freezing the slow variables.” We discuss the related problems 
later. The introduction of a small parameter /i formally ensures that the commutator [x, 7i] 
is small. There is no universal interpretation of this fact; this depends on each concrete 
physical problem. 

Let us also note that one can consider the equations from Examples 2 and 3 as infinite 
vector ones. To show this, let us assume, for simplicity, that, for each x G M”, the spectrum 
of the operator —{l/2)d‘^/dy‘^ + v{x,y) is discrete and simple and that the corresponding 
eigenfunctions {wn{x, y)} and eigenvalues A„(a;) depend smoothly on x. Then one can expand 
any solution \k(a;,//,t) of Eq. (I2.7jl in the Fourier series 


T = y^Wk{x,y)'ijjk{x,t). 

k 

Substituting solution dZHl) into Eq. (EH), we obtain: 


( 2 . 8 ) 


. dw„ di'n X 


2 dx'^ 




d‘^Wr. 

dx"^ 


fy'fn 


(2.9) 


If we introduce the inhnite vector = (//'i, ^> 2 , ■ ■ , then we can represent Eq. ()2.9j) as 

the following inhnite vector equation with inhnite-dimensional matrix Hamiltonian 7d: 

iyift = 'Hfi, id = Hq{p, x) + pHfip, x) + y^H 2 {p, x), 
f tP^ \ 1 cPiij 

{hdo)kn = f Y + An(a;) j {Ufikn = -i{wk, {hd2)kn = -^{wk, ~^)y 


In all examples considered above, the momentum operators corresponding to the slow 
variables Xj are —iyd/dxj. Of course, one can consider a general situation in which the 
Hamiltonian depends on the operators Xj, pj generating the Heisenberg algebra with com¬ 
mutators \pj,Xk] = ySj^k, /i <C 1. Such a situation appears in the electron-phonon inter¬ 
action, we shall discuss it in §4. The other obvious generalizations of the equations with 
operator-valued symbols are vector equations containing “slow” and “fast” variables. For 
instance, we can consider the Pauli equation in a thin quantum waveguide. In this case (see 
§4.6), the symbol is a matrix operator differential with respect to fast variables. 
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To conclude this section, we present a useful formula which plays an important role in 
the future consideration. Let A and B be pseudodifferential operators 

A = A(—i/ud/dx, X, ju), B = B{—ifid/dx,x,fi), 
then the symbol smb{AB) of their product AB is equal to (see j14p 

^ ^ 1 2 

smb{AB) = A{p —ifid/dx,x, ^)B{p,x, ^). (2-10) 


3 General scheme of the operator separation of vari¬ 
ables in adiabatic problems. 


3.1 General statement of the problem with operator-valued sym¬ 
bols and parameters. 

We are going to construct a certain asymptotic solution T = ..., s > 2, to vector 

equation (Q with a small parameter p, -C 1 or its stationary variant 

= ET. (3.1) 


In fl2.2j) and in ()3.1jl . the matrix operator (quantum matrix Hamiltonian) is generated by 
its operator-valued symbol 


Hii ... His 

n = n{p,x,-id/dy,y,t,IJ-) = . , 

Hsi ... Hss 


Hij = nij{p,x, -i—,y,t,y), I < i,j <s. 


(In the stationary case, Hij do not depend on time t.) We assume that the operator-valued 
symbol (the matrix-operator) H = \\Hij{p, x, —i-^,y, p)\\ smoothly depends on p,x,t and 
acts in an appropriate vector Hilbert space with coordinates y from some domain Mj^ and 
with the inner product (•, •)| 2 ^ (for instance, in L 2 (Mj^) x L 2 (My) x ... x L 2 (My)). Another 
natural assumption is that the symbol H{p,x, —i-^,y, p) can be expanded into a regular 
series with respect to the parameter p: 


d d d 

H{p, X, p) = Ho{p, X, A + yA) + ■■■■ (3.2) 


We also assume that the (pseudo)differential operator H acts in an appropriate expanded 
Hilbert space with coordinates {x, y) and all the future operations related to 

them are valid. Of course, one has to verify the last assumption in each concrete problem. 
Usually (but not always), we shall consider situations in which H, as well as the operator 
H, are essentially self-adjoint. 

It is important to emphasize again that in (El there is a small “adiabatic” parameter p 
in front of the derivatives with respect to “slow” variables x, but there is no small parameter 
in front of the derivatives with respect to “fast” variables y. 

Of course, one has to add additional boundary and initial conditions to Eq. ( 12 . 2|) . We shall 
do it later after the discussion in §4, and now we only note that we are going to consider only 
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special problems interesting from the physical viewpoint. The statements of these problems 
follow the adiabatic separation of the original Eq. (Q into a set of rednced eqnations 
corresponding to different “terms” or “modes” and determined by “effective Hamiltonians” 
or “dispersion relations.” We present onr concept of this separation (the “operator separation 
of variables”) together with the corresponding formulas in the two subsequent sections. 

3.2 Anzatz of the operator separation of variables. 

Let us illustrate the main ideas of the operator separation of variables in adiabatic problems 
with an example of a “quantum waveguide” (EH). If the potential v{x, y) is the sum ni(a;) + 
V 2 {y), one can separate the variables and hnd a special solution to Eq. (12.7|) as a product 
of two functions (modes) x(|/,/i)'^(a:, f, p). It is clear that this representation is not true if 
v{x,y) 7 ^ vi{x) + V 2 {y), nevertheless, since there are different scales in the longitudinal and 
transverse directions, we can separate the modes adiabatically. According to the standard 
adiabatic approach based on the fundamental papers by Born and Oppenheimer, the leading 
term of the wave function in the adiabatic approximation is sought in the form of the product 

T(a:, y, t, y) ^ y, f, y). (3.3) 

But this representation can be used in a situation when the function 'ip{x,t,y) is quite 
smooth and works poorly for large enough energies of longitudinal motion. If the func¬ 
tion '0(a;, f,/i) exhibits fast oscillations, for instance, if -0 is the WKB-solution '^(a;, f,/r) = 
exp{iS{x,t, y)/y)ip{x,t, y, h), then representation (I3.3|l is not convenient for the asymptotic 
expansion and, instead of formula (USD, one has to include the classical momentum dS/dx 
into the factor x(j/, a;,p.) and use the formula jH] 

T(a;, y, t, y) ^ x{dS/dx, x, y, y)ij{x, t, y). (3.4) 

Recall that the phase S is the solution of the Hamilton-Jacobi equation dS/dt -|- 
Hes{dS/dx,x,t) = 0 with the so-called effective Hamiltonian Hes{p,x,t). For the case 
in which Hes{p,x,t) is a function of p only and S = —uit + px, the Hamilton-Jacoby equa¬ 
tion is the dispersion relation. Formula (D is still not satisfactory, because for the case 
in which there exist focalization effects, i.e., there are turning points or caustics, the WKB- 
representation is not true, and it is necessary to change the form of V’ and the form of x- We 
propose to “correct ” (El) in such a way that a new formula would also work in the case of 
focal and turning points. This correction is based on the observation that, in the WKB-case 
modulo a small correction, the right-hand side in dH remains the same (see, e.g., [H]) if one 

1 2 

assumes that the hrst factor is the (pseudodifferential) operator x{—iyd/dx,x,y,t, y) writ¬ 
ten as a function (its symbol) of the noncommuting operators x and p = —iyd/dx. Finally, 
we suggest to look for the solution J/(a:, y, t) in the adiabatic approach in the following form 

mmi [131201: 


Q 2 


(3.5) 


where x is the “pseudodifferentional” operator whose symbol has the (asymptotic) expansion 
with respect to the parameter y 


x(p,x,y,t,ti) = x«(p,x,y,t) + pxAp.x.y.t) +.... 


( 3 . 6 ) 


From the physical viewpoint, representation (13-Sp means that we “freeze” not only slow 
variables x as in formula dH, but also slow momenta, which are differential operators 
—ijjid/dx in quantum mechanics. Note that in many situations the leading term Xoij>) Vi t) 
in expansion dSSl) does not depend on p, but the corrections usually do. This dependence 
plays an important role in the estimation of the limits of the adiabatic approximation in 
concrete problems. 

We still do not £x the equation for the function -0 describing the longitudinal motion. 
Following the idea of the so-called Peierls substitution in solid state physics (see, e.g., uu 
HH]), we assume that the wave function -0 is a solution of the following equation (describing 
the longitudinal dynamics): 


ip'ipt = L'lp, 


'' d 2 


( 3 . 7 ) 


where L is a pseudodifferential (sometimes, differential) operator with symbol L(p,x,t,p) 
having the expansion 


L{p, X, t, p) = Lo(p, X, t) + pLi{p, x,t) + .... (3.8) 

The operator L is called the (full) quantum effective Hamiltonian with the principal part Lq. 
Sometimes, the symbol Lq is also called an effective classical Hamiltonian and is denoted by 
Lq = Hes{p,x,t). The operator x will be called an intertwining operator (cf. [HI |2I1122) )• 
Equation dSIZl) can also be understood as the quantization of the Hamilton-Jacobi equation 
or the dispersion relation. The wave function has different names in different helds. 
For instance, is a. nuclear function in molecular physics, a longitudinal wave function in 
waveguides, an electron function in crystals, etc. We shall call it a wave function of adiabatic 
motion and we shall call Eq. dSIZl) the effective equation of adiabatic motion. 

Representation (13.5p together with Eq. (EZZl) (a generalization of the Peierls substitution) 
is a formalization of the operator separation of variables in the adiabatic approximation. Of 
course, the corrections Li, L 2 , ... appear in the problems in which the variables cannot be 
separated exactly. 

The reduced equation dSH) contains less independent variables and hence should be sim¬ 
pler than the original one. Thus, we see (and we mentioned this in Introduction) that solving 
the original equation can be divided into two parts: 1) the “operator {adiabatic) separation 
of variables” based on formula ()3.5|) . which reduces the original equation to Eq. dSUl), and 2) 
the process of solving this simpler equation. 

The realization of the hrst step consists in hnding the symbols (functions) Xj ci-nd Lj. We 
shall state the general scheme of their construction and discuss different related questions 
(e.g., concerning the reasonable number of terms in expansions (ESI), (EHl)) in the next 
sections. Now we discuss a natural generalization of the operators x ^-nd L. 

It is easy to see that, in the case of exact separation of variables, y = xo is an eigenfunc¬ 
tion of some additional spectral problem. The same fact holds for the functions Xo{x,p, y, t); 
later we shall numerate them by a multiindex u. Thus formula (13.5|) describes only some 
special solutions of the original equation corresponding to the term with the index u. It 
is possible to construct more general ones summing solutions (13. 5p with different indices u 
and the corresponding y, fj. Another conclusion is that, in the case of exact separation of 
variables, the spectrum of the above-mentioned additional spectral problem can be degener¬ 
ate and several eigenfunctions can correspond to the same eigenvalue. Then, instead of the 
product y(|/,/i)'^(a;, f,/i), one should write the sum Xi(2/, h)V’i(2^, h), where k is the 

multiplicity of the corresponding eigenvalue. The same generalization should be performed 
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in formula (ESI). Also if the original problem is a vector one (i.e., if (EH is a system of 
PDE for s unknown functions), then Xj has s components. Finally, in formula and in 
Eq. ()d.7j) we mean the following: 

1) X is a matrix pseudodifferential operator with s rows and r columns, 

2) Ip is an fc-dimensional vector function pj = {'ipi,..., 'ipkY', 

3) L is an fc X A: matrix pseudodifferential operator with the principal symbol Lq = 
Hes{p,x)Ek, where the number r determines the multiplicity of the corresponding effective 
Hamiltonian iJeS and Ek is the k x k identity matrix. The corrections Lj usually are not 
diagonal, which means that interaction is present inside the mode (or the term) determined 
by this effective Hamiltonian Hes{p,x). 

The number of terms in the expansions of the intertwining operator x and the operator L 
(with hxed index z/) can be arbitrarily large. However, it is, as a rule, a very complicated 
problem to calculate the terms of these series explicitly, even terms with small numbers. 
Therefore, it is natural to consider only the terms required to estimate correctly the leading 
term of the asymptotics of the wave function or of the energy value. 

However, the notion of the “leading” term of an asymptotics can be determined not 
only by the adiabatic parameter p, but also by the other ones, for instance, by the so- 
called “semiclassical parameter” h, which is related to the form of the effective potential and 
the solution of the reduced equation (EH). The appearance of this new parameter is very 
important for future constructions of the asymptotics. We shall discuss the corresponding 
questions in detail later in §5. Now we only say that, usually, for the construction of the 
leading term of asymptotic solution, it is sufficient to hnd Lq, Li, and L 2 |p=o- Another 
interesting fact is that the effects of a semiclassical splitting of the effective Hamiltonian and 
a change in the classical characteristics occur in the degenerate case (see §4.6). 


3.3 Scheme of the operator separation of variables 

To simplify the future consideration, let us assume that, in Eqs. E3) and EH), H and the 
operator are essentially self-adjoint. We shall seek the solution of Eq. EH in the following 
form: 


'^i{x,y,t,p) = '^Xijik 


d 


(3.9) 


where pj = {'ipi ,..., is the wave function of some chosen term (or a chosen “fast” mode) 
with the degeneration multiplicity equal to k and x is an intertwining matrix pseudodiffer¬ 
ential operator: 


Xii • 

• Xlk 

Xsl • 

Xsk 


X(p, X, y, t, p) = xo(p, X, y, t) + pxiij), x,y,t) + .... (3.10) 


We assume that the vector function -0 satishes the “effective equation of adiabatic mo¬ 
tion” ()3.7j) generated by the matrix operator L 



L{p, X, t, p) = Lo{p, X, t) + pLi{p, x,t) + ..., 


where the matrix Lo(p,x,f) is proportional to the unitary k x k matrix Ep. Lo(p,a;,t) = 
HeffEk. The coefficient of proportionality iLeS is an effective Hamiltonian. Hence the problem 
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is reduced to finding the operators x and L or their symbols x and L. After we find them, 
we can reduce the initial problem to a more simple (reduced) equation ()d.7j) for the vector 
function ip. The original solution T can be reconstructed in accordance with dSSD. 
Substituting the function T from (USD into Eq. (E3), we obtain: 

Using condition dSZZD, rewrite this equation in the following form: {xL + iuxt — 'Hx)'^ = 0. A 
sufficient condition for the last equality to be valid is the operator relation xL+ifiXt—Hx = 0- 
Let us pass from operators to symbols [0] in this relation using formula (ITTUl) . This leads 
to the equation 


0 2 

x{p n)L{x,p,t, n) + inxt{p,x,y,t, IJ,)- 

1 

d 2 d 

-n(p ix)x(p, X, y, t, p) = 0. 


(3.11) 


It can be solved using regular perturbation theory, i.e., expanding the items into series 
with respect to p. Collecting terms of order fjP = 1, we obtain a family of spectral problems 
for the self-adjoint operator T-Co{p,x,y, —id/dy,t) depending on x,p,t: 

d 

Uoip, X, -f—, y, t)xo{p, X, y, t) = xo{p, x, y, t)Lo{p, x, t). (3.12) 

dy 

We shall assume that the asymptotics (USD is completely determined by the eigenvalue (term) 
Hes{p,x,t) whose multiplicity k does not depend on p,x,t. Moreover, we shall assume that 
the value Hes is separated from the other eigenvalues or a part of the spectrum of T-Lq (if the 
spectrum contains a continuous component) uniformly with respect to {p, x, t) in a certain 
fixed domain {p,x,t) G A4. 

So 

Lo(p, X, t) = Hes{p, X, t)E, (3.13) 

where E is a unitary k x k matrix. The matrix Xo{x,p, y, t) consisting of orthonormal vector 
columns, i.e., eigenfunctions of the operator Hq corresponding to the eigenvalue ifefr(x,p, t), 
is the intertwining operator on the proper subspace induced by this eigenvalue. It is natural 
to assume that Xo{x,p,y,t) depends smoothly on all its arguments. 

Collecting terms of order /i, we obtain inhomogeneous equations for Xj and Lj: 

(Tio - HesE)xj = Ej - HjXo + XoLj, j = 1,2,... (3.14) 


where Ej depend on xo, ■ ■ ■, Xj-i and Lq ..., Lj_i, in particular, 

' dUi dxo _ dxo dLi 
dpj dx^ dpj dx^ 


Fi = Pxo, F2 = Vx\ - niXi + XiLi + I Y, 


+ 


Here 


V = i 


A 

'Wt 




8‘x'i d^xs 




dpidpj dx^dx^ dx^dxl dpidp 


3 2 




dHo d d 


dpj dxi dxi dpjj 

_ _ 
dt dt 


= i 


A 

dt 




dHo OH, 


eff 


d ^ d 

^ dx^ dpj 


E 


3 

d 

dpj dx^ 


dpj dpj 


d 

dx^ 


(3.15) 
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Due to the self-adjointness of the operator {TiQ — He^E) and the Fredholm alternative, the 
solvability condition for this equation is equivalent to the condition that its right-hand part 
is orthogonal to the vector-columns of the matrix xo- If follows that Lj = ,'HjXQ)\y — 
(XQ,Fj)|y. In particular, one can obtain: 


Li = 




dxo \ 

dt / 



dUo 

dpj 


dHes ' dxo 
dpj dx^ 


y 


(3.16) 


Assuming that Li has form (jd.lbjl . one can find the correction, i.e., the matrix Xi = 
{Ho — HesE)~^{Fi — HiXo + Xo-hi), fixing it for determinacy by means of the condition 
of orthogonality of vector-columns of the matrices xo and xi- Th® reiteration of this pro¬ 
cedure leads to calculation of Lj,Xj- Formulas (Id.ldj) . (Id.lbj) . etc. give the coefficients of 
expansion of the symbol of the reduced equation (Ell). Note that the construction of the 
correction Li includes only functions of zero approximation (like in the standard perturba¬ 
tion theory). In general, the symbol L 2 includes Xi, so, to find it, one has to invert the 
operator {Hq — H^sE). 

Remark 1. The methods from ^1] allow one to consider more general situations in 
which the quantum Hamiltonian H{p, x, Py, y, p) is a function of vector operators {p, x, Py, y) 
with commutator relations [%,Pj] = ip, [yj,Pyj] = i, /i -C 1, or even more complicated ones 
(see Example in Section 4.2). However in this paper we basically consider the situation in 
which X = x,px = —ipd/dx. 

Remark 2. It is not difficult to modify the presented formal scheme for the non¬ 
self-adjoint original operator H. In particular, one has to use the eigenfunctions of the 
adjoint operator in the orthogonality conditions. But, of course, it is necessary to add 
some additional conditions like the existence of real-valued effective Hamiltonians, etc. (see 
Example in section 4.3). 

Remark 3. Operator separation of variables and adiabatic approximation in 
classical mechanics. There exist a certain classical analog of adiabatic approach based on 
the “operator separation of variables” (1231211123123 EH). The main idea can be illustrated 
by means of the Hamiltonian H{ppx,x,py,y, p) with a small parameter p. If we change 
the variables x, by ^ = x/p and p^ = ppx, then we obtain a Hamiltonian of the form 
H{p, p^,Py,y, p). It is convenient to write the Hamiltonian in noncanonical variables x and 
p, dp A dx = pdpx A dx = pdp A d^: H{p,x,py,y, p). The Hamiltonian equations for the 
variables p, x, Py, y have the form 


on , 

X = p— < 1, 


on , 

p = < 1 , 

ox 


y = 


m 

dpy 


Py = 


m 

dy ■ 


(3.17) 


Since we have p, x ^ p for the derivatives, while p^, p ~ 1, it is natural to say that the 
variables p, x are “slow variables” and the variables Py, y are “fast variables.” Taking into 
account that there are variables of two types, it is natural to “freeze” slow variables and 
obtain a family of Hamiltonians with k degrees of freedom depending on the parameters 
(p, x). We do not consider the resonanace problems and restrict our consideration to the 
case in which k = 1 and {py,y) G Let us assume that, in some domain (p, x) G D, the 
trajectories of H{{p, x},py, y, 0) are closed. The braces {■} mean that the included variables 
are considered as parameters (i.e., are “frozen”). Then it is possible to introduce “action- 
angle” variables (J, <p) corresponding to these closed trajectories. The passage to these 
variables is determined by the change of variables y = Yo{J,(p,p,x), Py = Py{J,ip,p,x). 
Unfortunately, this change of variables is not canonical and, to make it canonical, one has 
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to add corrections and write y = lo(J, ip, P, X) + ip, P, X) + ..., py = Py{J, p, P, X) + 

pPy{J,p,P,X) + ..., p = P + pPi{J,p,P,X) + ..., and x = X + pXi{J,p,P,X) + .... 
Then the original Hamiltonian can be written in the form 


n{j,p,p,x) = noUP,x) + p 


[P, X, F„", Y„)y\ + ^(F, A', F», Y„)P^+ 


dy 


dpy 


+ 


dUo 


dn, 


{P,X, yo)Xi + V^(P,X,P° Yo)p^ + n,{p, X, p° Xo) 


+ 0(/i^), 


P^ = P^{J,p,P,X), Y, = Y,{J,p,P,X), P^ = P\j,p,P,X), Xi = Xi(J,v?,P,X) 

P; = P°(J,¥.,P,X), Xo = Yo{J,p,P,X), (3.18) 


where T-[{p,x,Py,y, p) = T-Co{p,x,Py,y) + pHi{p,x,Py,y) + 0{p^). Now we have a typical 
problem from averaging theory. After averaging, we obtain two terms of the expansion of 
the effective Hamiltonian L(J, P, X, p) = Lo(J, P, X) + pLi{J, P, X) + 0{p‘^): 


LoUP,x) = noUP,x), 


p2'K 

Li(J,P,X)=/ , 

Jo 


dp 


>0 

dUn,^ ..... OK 


X, F», Y,)Y\ + ^(P, A, P», F„)F>+ 


(PA'.fJ.yy.Yi + -^(p.Y.pJ.yyP + Hi(pA',p“,y„) 


(3.19) 


so J = const and the integration of the original system is rednced to solving a system with 
n degrees of freedom. The action J corresponds to the “qnantnm” nnmber v of the term x'^i 

the term dpTLi ^P, X, Py{J, p, P, X), Xo(J, p, P, X) j corresponds to (xo, 'HiXo)yi and the 
corrections related to the canonical change of variables in the classical problem correspond to 
other terms in Li in qnantnm problem. Of conrse, this is simply an analogy (cf. [71 12811^ . 
etc.). 

Remark 4. The classical analogne of the redaction at the first stage is well known 
EZ|: exclnding the fast variables, we obtain a system in the zeroth approximation with 
holonomic constraints; this system is eqnivalent to the n-dimensional Lagrangian system. 
Thns the adiabatic redaction to Eq. (Ell) conld be interpreted as the “exclnding of qnantnm 
constraints” (see 1201 ED E2 ESI El E3 Ei ) • Bat the classical system corresponding to the 
rednced qnantnm system generally does not coincide with the resnlt of the classical redaction 
(in the sense EH). Moreover, the classical systems arising in the adiabatic redaction tarn 
ont to be different for different relations between p and “semiclassical” parameter h, which 
will be introdnced in §5; nsing the terminology in in EH], we can say that they correspond 
to different nonstandard characteristics of the qnantnm problem. For example, the classical 
eqnation of motion in nanotnbes can sometimes inclnde terms arising becanse spin exists. 
This is explainable physically because longitudinal motion is already determined by rather 
small energies that are quite comparable to the spin energy. 

Remark 5. In some problems one can apply the semiclassical approximation to solve 
Ea. ()3.12j) (see e.g. iniEnj). 


4 Examples of problems with operator-valued symbols 
and parameters. 

Let us illustrate the general scheme with several examples. 
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4.1 Equation with rapidly oscillating coefficients and electron 
waves in crystals 

One can meet equations with rapidly oscillating coefficients in many problems of solid physics 
and continuum mechanics. For instance, these types of equations describe propagation of 
electron waves in crystals, of elastic waves in composite materials, etc. For constructing 
asymptotic solutions to equations with rapidly oscillating coefficients, there exist different 
approachs adapted to their behavior and their properties. Among these approaches, one can 
at least mention averaging methods, homogenization, and adiabatic approximation. There 
exists a very extensive literature concerning this topic and a review of all these approaches 
is very far from our aims. We only want to show that one can look at the equations with 
rapidly oscillating coefficients from the viewpoint of equations with operator-valued symbols 
nnna and, for the construction of their asymptotics, use the method described in §4. So 
here we mention only the general monographs P El EOl il EH 112 and ideologically close 
papers US El EH- 

As an example, we consider the Schrodinger equation with fast oscillating potential 

ijji'ipt = , x)'ip, X e M”, $ e R”^, (4.1) 

^ jJj 

where v{y,x) is a smooth function 27r-periodic with respect to each “fast” variable yj, j = 
1,... ,m. The given phases 4’j(x) are smooth functions. Generally speaking, their number 
k can be arbitrary. In some problems, the phases are linear functions = {kj,Xj)] the 
case of nonlinear phases describes the case of a nonuniform potential v. The additional 
dependence of the potential on the variable x implies its slow deformation. Eq. (EH) with 
such a potential simulates the propagation of electron waves in a lattice or, for instance, if 
m = 1, in stratified media. The simplest example of a rapidly oscillating potential is given 
by the formula (n = l,m = 1): v = Vo(x) -|- a(x) cos where Vo(x),a(x) are smooth 
functions. 

Let us find the unknown function ip{x,t,y,) in the form 

V’(x, t, /i) = d' X, t, , (4.2) 

where the new unknown function 4/(?/, x, t, y) is 27r-periodic with respect to each variable yj. 
Substituting EH into Eq. EH, we see that the function 'ip{x,t,y) ()4.2|1 satisfies Eq. EH 
if the function '^{y, x, t, y) is a solution of Eq. (I2.2|l with 

/ d d\‘^ 


We again see that the small parameter y is in front of the derivative d/dx, but there is no 
small parameter in front of the derivative d/dy. Thus the equation with Hamiltonian (14.dj) 
is an equation with operator-valued symbol, namely. 


n = no- 

i=i 


^dy 


T<o= P - * 


j 




i=i 


d 


dx dyj 


+ x). 


(4.4) 


For each fixed (p, x), the operator Ti acts in the L 2 -space on a fc-dimensional torus. 

If the vectors are linearly independent for each x, then the spectrum of the operator 
7i is discrete, but, to obtain the reduced equation, one must also use the assumption that the 
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multiplicity of eigenvalues is independent of {p,x). In particular, if m = 1 and d^/dx 7 ^ 0, 
then the spectral problem (IXT^ for determining the effective Hamiltonians is a periodic 
problem and can be reduced to the 1-D Schrodinger equation on a circle for Bloch solutions 
(see iiniiiEi). To realize this reduction, we change the variables in the equation T^oXo = 
HesXo Xo cind -ffeff as follows: 


y = Ui, Xo = ue 



and put 


P = P{p, x) 


Then this equation takes the form 



(4.5) 


—U^^ + v{U^,x)u = Eu, S = PleS — p^ + P"^, 


and the periodicity condition becomes the Bloch condition: 

x) = x). 


The variable x is contained in the reduced problem as a parameter. The variable (number) 
P is called the quasimomentum of the corresponding Bloch solution. It is a well-known fact 
that the spectrum of the operator —d‘^/d^‘^ + v{U^^ x) on a circle consists of bands and gaps. 
Let us enumerate the bands by the number v and denote the ends of the z/th band by and 
The spectral parameter £ and the quasimomentum P in each z/th band are connected 
by the dispersion relation 

£ = £'^{P,x). 

The assumption on the potential v{y,x) that the z/th effective Hamiltonian (eigenvalue) of 
the operator PC^ is simple (or does not intersect with other effective Hamiltonians at some 
points {p,x)) is equivalent to the assumption that, for each x, the uth band does not stick 
together with the z/ — 1st and z/-|- 1st bands. Under this assumption, one can hnd the reduced 
equation describing solutions corresponding to the z/th term (effective Hamiltonian) 

x),x)+p^- (P(p, a;))2. 


The corresponding function XoiViVi x) is expressed via the Bloch function P, x), by 

the formula 


5 4 / P{p,x)y. 

Xo = u {—,P{p,x),x)exp{-i - - —, 


If n = 1 and <I> = x, then P = p, which leads to the well-known fact in solid state physics: 
the quasimomentum becomes the momentum for the equation for electron waves in crystals. 
In contrast to Examples 1, 2, and 3, the effective Hamiltonian here is not a polynomial in p 
and the function Xo depends on the momentum p. 

Let us write the hrst correction Li to the effective Hamiltonian. Using formula we 

obtain 


Li = (xo, + (Xo, (2p - -^)^Xo) + 2iRe(Vxo-J^, Vxo)- 


Note that it is possible to meet a situation in which the number of phases in the potential 
is greater than the dimension of the conhguration space n. For instance, consider the case 
m = 2, n = 1. In this situation, the operator PL is degenerate and its spectrum has a 
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rather complicated structure. In particular, the spectrum can be an everywhere dense set 
on the spectral axis (the so-called devil’s stair) and the multiplicity of its eigenvalues can 
depend on x. These problems are related to different types of complicated resonances and 
to problems with intersecting characteristics. Some results for the solutions of the Cauchy 
problem in this situation are obtained in Ha 1201 Hi. 

Remark. Bloch electrons in a weak magnetic field and the Peierls substitution. 
The original Peierls substitution was hrst proposed for the problem with Bloch electrons in a 
weak magnetic held (see, e.g., j3]). Using the terminology of this paper, this problem can be 
stated as a spectral problem for the magnetic Schrodinger operator with the periodic electric 
potential: 

+v(^}, (V,A>=0. 

For simplicity, we restrict ourselves by simple cubic lattice, i.e., assume that v = n(|/i, 2/2,1/3) 
is 27r-periodic with respect to each variable yj = Xj/y, j = 1, 2, 3. The fact that the magnetic 
held is weak means that A{x) does not contain any irregular dependence on the parameter y 
in contrast to the crystalline potential v(^). After the space regularization dOl) with phase 
vector <F(x) = x similarly to ()4.3j] . we obtain the problem with operator-valued symbol 
7i = (1/2)(p — id/dy — A)^ -|- v{y). The procedure based on the formulas (I4.5|l leads to the 
problem for Bloch solutions for the operator Ti = —(l/2)Aj^ -f- v{y). Let £ = £’^{P), P = 
{Pi, P 2 , P 3 ) determine the dispersion relation for the Bloch solutions with quasimomentum 
P. Then, according to this subsection, Lg = = £''{p — A), U{ = [tr(^^|^)](p — A), 

and in (j3.7j) L = Lq + yL^ + 0{y?) = £’^{—iy-^ — A) -\- 0{y£‘) which is exactly the Peierls 
substitution. The higher order corrections look much more complicated than the leading term 
BH. It seems to us that our approach allows one to calculate these corrections easier than 
in BH. Semiclassical analysis of the reduced equation (EZD reveals a very complicated and 
striking topology of surfaces invariant to the corresponding phase flow. The recent results 
and bibliography can be found in [IH] . 

4.2 Electron-phonon interaction 

As was noted in Remark 1 in Section 3.3, §3, one can consider adiabatic problems as problems 
containing “slightly noncommuting” operators. In the zeroth-order approximation, these 
operators can be substituted by “c-numbers,” which allows one to determine a term. For 
“slightly noncommuting” operators, there are physical quantities slowly varying in time. The 
slightly commuting operators can generate a certain Lie algebra. In the simplest case, this Lie 
algebra is the Heisenberg algebra, and we can directly use the scheme and formulas proposed 
in §3. Sometimes, it is possible to consider the same problem from a different viewpoint, 
which depends on the choice of the operators. For instance, to simplify the form of the 
original Hamiltonian, one can use a (noncanonical) change of variables, which, in turn, leads 
to the replacement of the Heisenberg commutation relations by different ones. We consider 
such an example which can be analyzed from these different viewpoints, but since the study of 
problems based on non-Heisenberg commutation relations requires new nontrivial algebraic 
and geometric constructions, we restrict ourselves to the approach described in §3, although, 
from some viewpoint, the approach based on non-Heisenberg commutation relations can be 
more readily realized in some concrete problems. 

The electron-phonon interaction is the interaction between light fermions (electrons) 
and heavy bosons (phonons). Here the lattice modes (bosons) are slow and the electrons 
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(fermions) are fast. The Hamiltonian of electron-phonon interaction is (see, e.g., |49p 

'H = + i^n'^n+l) + ^K{Xn+l - Xnf + pI/2, (4.6) 

n n 

[Xn,Pn'] = ipSnn', P = H/ay/moCQ (4.7) 

'0fc]+ = '4’t'^k + = iSjk- 

Here mo, cq, K, a are physical constants; moreover, /i -C 1; the linear operators ..., 
and Ti,..., act on the Hilbert space Hi, the linear operators Pi,... ,Pm and hi,..., xm 
act on the Hilbert space H 2 . The full quantum Hamiltonian acts on the Hilbert space 
H = Hi 0 H 2 . The typical situation is given by the operators 

Cn = f{Xn+l - Xn), (4.8) 

where f{z) is a smooth function; in particular, = 1 — a(h„_|_i — Xn)- Let us set X = 
(hi,... ,hAr), P = (pi, .^.. ,Pjv), T = (-01,.. . ,-0^1), In a more general 

case, the Hamiltonian ?h of electron-phonon interaction can be written as 

7h = (T+,/:(X,P)T)+ <h(X)+ PV2, (4.9) 

where C{P,X) is a Hermitian M x M matrix with coefficients depending on X and P, and 
the operator C{X,P) is understood in the sense of Weyl calculus (see [HI [HI I221)- To 
simplify the consideration, we restrict ourselves to the case in which C, = CpX) -|- P 2 (-P); 
then the question about the ordering of the operators X and P does not appear. As we 
have just said, it is possible to develop the “operator separation of variables” based on the 
algebra of operators Cj,pj with the commutation relations 

[Pj, Cfc] 5jk+pCk, 

but here we use the standard representation Xj = Xj and pj = —ipd/dxj and close the lattice 
by the Born-Karman periodicity condition identifying the operators with the numbers j and 
j + M. 

Denote the electron-phonon wave function by T and consider the stationary problem 

HT = ST, T G H. (4.10) 

We have an equation with operator-valued symbol which, obviously, is the operator 

n{X, P) = (T+, C{X, P)T) + ($(A) + PV2)/i (4.11) 

acing on the Hilbert space Hi. We also denote the identity operator acting on Hj by Ij. 

To realize the scheme of the operator separation of variables, it is necessary to hnd 
the spectrum of the operator-valued symbol 'H{X,P). A nice fact is that this spec¬ 
trum can be expressed via the eigenvalues of the matrix C. Namely, suppose that 
... ,'ipMiEj)) are the eigenvectors of the matrix C{X,P) corresponding 
to its eigenvalues Ej {X, P) Ei < ... < Em and satisfying the normalization conditions 
{'ip{Ej),'ip{Ek)) = 5jk- Using the basis {'^(H^)}, one can expand the operators T+ and T 

M M 

'il!+ = = '^ij{Ej)dj. (4.12) 

i=i j=i 
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The coefficients (operators) of this expansion, and dj, are called the creation and annihi¬ 
lation operators They determine the operators of the number of particles Afj = d'^dj. 
Finally, the operator-valued symbol can be written as 

M 

n{X, P) = Y, + ^V2, = apj. 

i=i 

From this, we hnd the z/th effective Hamiltonian 


P) = Y, P) + 4(A') + PV2, (4.13) 

i=i 

and the symbols of the corresponding intertwining operators Xo'- 

= ATi ... AV) ■ 1. (4.14) 


Let us analyze the formula for Li. If the effective Hamiltonian ifgg is degenerate or 
Ci{X),C 2 {P) are complex matrices, then one must use the general formula (ld.lfi|) . For 
the case in which C 2 {P) = 0 and C.i{X) is a real-valued matrix and its spectrum 
(i?i(x),..., Em{x)) is nondegenerate, it follows from formulas ()4.1d|l that the effective Hamil¬ 
tonians are nondegenerate and and Xo Thus, taking into account the relations 


^1 = 0 , {xI,^)y = ^Tt{Xo,Xo)y = 0, ^ 


dHeff J _ 


75 —/ = 0, we obtain 

dpi ’ 

Li = 0 . 


This equality holds for the operators c„ given by formula (gSl), in particular, if c„ = e*" 

It was pointed out that in this case it is possible to relate problem (ICTIl) to the integrable 
Toda lattice model IsilEg. Then this observation was used to construct its semiclassical 
asymptotics insi, and at the same time, commutation relations of different type were chosen. 


4.3 Internal waves in ocean in a pycnocline 

The next example concerns a situation in which the operator Ti is not self-adjoint, moreover, 
the system of equations under study differs a little from (Q. Thus the scheme of §3 requires 
a slight modification. 

We consider a system of hydrodynamic equations for an ideal incompressible liquid lin¬ 
earized on the flow with velocity U and density pq. Let x = {xi,X 2 ) be the horizontal 
coordinates, z be the vertical coordinate, u = (mi,M 2 ,M 3 ), and p be perturbations of the ve¬ 
locity and density, H be the pressure, g = (0, 0, | 5 f|) be the gravity acceleration. We introduce 
dimensionless variables and parameters by the formulas: U = f/'wiA, u = u'oJiX^ t = t'uj 2 , 
X = x'L, z = z'L, g = g'ujfX, H = n'a;f A^, and po = pgp, p = p'p, where oji is the character¬ 
istic frequency of the internal wave, A is the characteristic wavelength, L is the characteristic 
distance in the horizontal direction within which the characteristics of the liquid vary, 1 x 2 is 
the average value of the Vaisala-Brunt frequency, p and are the characteristic values 
of density and pressure, and p = XjL = uji/ui 2 ^ 1 is a small parameter. 

In the dimensionless variables, the linearized system for waves in liquids has the form 

HOI El: 

f + l^Po{U, V)m -h PPo(m, V)f/ -h pVH -H pp = 0, 

I pf+ p(17,V)p + /i(M,V)po = 0, (4.15) 

I p(V,m) = 0, 
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where {•, ■) is the inner prodnct in The difference between Eqs. and Eq. (Q is 

that the last eqnation of this system contains the time-derivative. 

We shall consider that the width of the picnocline Az ~ p, and that the depth of its 
location varies at distances ~ 1, i.e., po = n — f{x),x) (the equation for the surface of 
a picnocline is z = p/). The functions U{x) = {Ui,U2,0), po{y,x), f{x) are assumed to be 
smooth, 0 < ^0 < Po < Co (^O) Cq are constants), the square of the Vaisala-Brunt frequency 
^0 — “IpI^/Po is positive and vanishes rather fast as \y\ —>• oo. We assume that the 
functions Uj and p decay quite fast as \z\ —> oo. The other boundary and initial conditions 
for Eqs. (g:T3ii are chosen in a special way and should be formulated for the corresponding 
reduced effective equation of adiabatic motion (uni). 

Let us introduce a new independent variable y = z/p — f{x) and a vector with Eve 
components T(a;, y, t, p) = (n, p, 11). Into the equations for the vector T, we must substitute 

the differential operators as —*• -and p-^ —>■ -S-. Then, for T, we obtain a 

system of equations containing “fast” variable y and slow variables x, t. For convenience, 
we multiply this system by z = \/ —1: 


ipB'^t 




. d 


p)T(a;,z/,t,p), 


B 


/po 0 0 0 0\ 

0 po 0 0 0 

0 0 Po 0 0 

0 0 0 1 0 

\0 0 0 0 0 / 


(4.16) 


The matrix symbols Tio and Tii have the form 


= -*Po 


T-Co{x,p,y, 

/8U1 

8xx 

8U2 

8U1 

8x2 

dlJ2 

8x1 

dx2 

0 

0 

0 

0 

V 0 

0 



0 0 0 \ 
0 0 0 
0 0 0 
0 0 0 
000 / 


[po{U,p) 0 0 0 Pi \ 

0 po{U,p) 0 0 p2 

0 0 Po{U,p) -i\g\ -i^ 

0 0 -i^ {U,p) 0 

V Pi P2 0 0 / 

/ 0 0 0 0 


+ 


0 

0 

jdpo I 7 df dpo 
''8x1 ^ ' 8x1 dy 

\ M-M) 


0 0 

0 0 


: dpo I -^dpo n 
' 8x2 ' ' 8x2 dy 

0 


0 

0 

0 

0 




M. 

8 x 1 

M. 

8x2 



0 

0 

0 


+ 


/ 


+(CV/)s(>|). 


where p = (pi,P2,0). Let us show that the general scheme of the operator separation can 
be easily modified for this situation, although one of the equations in the system does not 
contain the time-derivative. 

We seek the solution in the form (Hi) and assume that 'ip{x,t,p) satisfies an effective 
equation of the form (Q. Then, instead of (IXTTl) . we obtain the following relation for 
symbols: 


2 d 

Bx{.x,p -ip—,y,t,p)L{x,p,t,p) + ipBxt{.x,p,t,y,p)- 


2 d d 

n{x,p -ip—,y,-i—,t,p)x{x,p,y,t,p) = 0. 
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Expanding this relation into a series with respect to a small parameter /i, instead of the 
eigenvalne problem, we obtain the problem for the spectral parameter H^ff of the operator 
sheaf {Hq - BHes), 


'Hoix, p, y, t) - BHes{x, p, t)^ xoix, p, y, t) = 0, 

and the following eqnation for the corrections xi ci-nd Lp. 
d \ 

'Ho{x,p,y,-t—,t)-BHea{x,p,t)]xi = Fi-niXo + BxoLi, Fi = Vxo (4.17) 


d 




dHo d 


B- 


-dH^a d 


dpj dx^ dx^ dp 


.^d , .^diUo-BH,^) d 


dpj 


dx^ 


We assnme that the chosen spectral parameter Hes{x,p,t) is nondegenerate, so xo is a 
vector with Eve components {xoj}, j = 1, • ■ ■, 5. Then the solution of the problem for the 
operator sheaf has the form 


_ -pi dw{y,x,a) 


XOl Qy ■) 

.. ^P2 dw{y,x,a) 

X02 p2 Qy 1 

X 03 = W {y,x,a ), 

, , _ ,• dpQ w{y,x,a) 

T04 ^ Qy ^ , 

A dw{y,x,a) 


(4.18) 


X05 — 


dy 


Here A = {U,p) — Hef[{x,p), a = A^/p^, p^ = p\ + pi-, w{y,x,a) is an eigenfunction of the 
problem 

\ d d iy x^ 

-7^Po(l/, x)—w{y, X, a) + ° ’ w{y, x, a) = k{x, a)w{y, x, a), 

Po{y,x)dy dy a 

and Hefi{x,p,t) is a solution of the algebraic equation 

f {{U,p) - Hesix,p)f\ 2 
- )=P- 

We choose an eigenvalue x and consider the corresponding function Hes{x,p). In general, 
this function is multi-valued; we £x one of its branches and assume that this branch is a 
smooth function of x, p. 

The first correction Li is found from the solvability condition for Eq. (gUD. Its right- 
hand part must be orthogonal to the kernel of the adjoint operator {TL — BHes)*- We denote 
a function from its kernel by Xo- Then the hrst correction to the effective Hamiltonian is 


Li{x,p,t) = 


iXo: Bxo) 




n-/ d{Ho — BHes) d 


dp, 


dx^ 


Xo 


To determine xh we note that the construction of the operator adjoint to (Hq — BHes) 
is equivalent to the replacement \g\ —dpo/dy. This gives 


X*ok = Xok, /c = l,2,3,5, xm = -'i\9\ 


w{y,x,a) 

A 


20 






















Using these relations, we obtain 


(Xo) Bxo)y — Po 




dw 


dy 


dy+ po\w\ dy + 


PO^O |„..|2 


dUi 


{x*o,'HiXo)y = 

j,k=l 


P" 


-4(pW/> 


dw d f dw\ , \ 


+^(f^.V/) / P« 




dy dy 
dw dPw 
dy dy'^ 


PjPk I Po 
1 


dw 


dy 


A2 


dy- 


w\ dy, 


(4.19) 

(4.20) 


_ dw 


^rdpo \ _dpodw 
p, Vpo - V/— ) w——dy+ 


\ 9 \ 


dy / dy dy 
d f dpo 


dy + i{U, V/) / pow—dy - i{U, V/)— / w— —w dy. 


dy 


A2 y dy \ dy 


Using formulas (juisii-doni) we can calculate the hrst correction Li. We don’t give here the 
explicit formula for Li in the general case because of its bulk. 


4.4 Electromagnetic waveguides, integral optics, surface gravity 
water waves and shells. 

The electromagnetic wave propagation in waveguides is described by wave equation con¬ 
taining the second time-derivative. In the two-dimensional case, we have a situation similar 
to that considered in Example 2, §2. One can easily generalize the scheme of §3 to this 
situation. The change consists in the following: instead of one must write the second 
time-derivative in the left-hand side of (I3.7|l . The same change allows one to con¬ 

sider the three-dimensional waveguide problems. But now it is possible to consider waves in 
thin films (integral optics) or in thin tubes. The stationary variant of such equations is the 
Helmholtz equation 


(A-I-/c^n(a;))4/= 0 (4.21) 

with the refractive index n{x) and, e.g., the Dirichlet conditions T = 0 on the boundary of 
the film or the tube. The parameter p characterizes the ratio between the transverse and 
longitudinal dimensions of the waveguide; one can apply the adiabatic approximation if the 
boundary of the waveguide changes slowly. V. P. Maslov considered problems of this type 
in 1958 in [^, where he constructed asymptotic solutions and predicted the possibility of 
the construction of one-mode resonators by means of the waveguide geometry. Later on, 
problems of this type were considered in more general situations in optics and quantum 
mechanics (see, e.g., ISH E2l ESI El E3 Eni ESI EZl Ei Ei ) • 

The consideration of the Helmholtz equation is very similar to the consideration of the 
stationary Schrodinger equation. In what follows, we consider the Schrodinger equation in 
a quantum waveguide in a more complicated situation. 

More complicated examples similar to planar waveguides (thin films) lead to problems 
about the wave propagations in shells. Instead of the wave equations, one must consider the 
Lame equations in elastic theory. The operator separation of variables can also be used in 
such problems, but the study of these problems is far beyond the aim of this paper. 

A more exotic example of the operator separation of variables is given by the theory of 
surface gravity water waves over an uneven bottom (see, e.g., ISDl). Actually, this problem 
is the linearization of the problem with free boundary and the anzatz (13.5p was hrst used in 
this situation HH. The operator approach is discussed in detail in HH12111221 EH, so we do 
not consider this problem here. 
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4.5 Nanophysics: wave propagation in nanofilms. 

The two subsequent examples (quantum waves in nanofilms and nanotubes) are probably 
the simplest ones in adiabatic problems. An interest in these problems appeared recently 
because of great progress in nanotechnologies. It seems that most of the results described 
below, as well as many recent mathematical results (e.g., isiiis2iEsiisniinziEHiini), did not 
appear many years ago, because there was no deep physical interest in the corresponding 
problems. Now the question is to study concrete applied problems. Needless to say that the 
representation of the solution in the form appropriate for practical analysis is an additional 
and sometimes nontrivial problem (even in the case of nanofilms and nanotubes with simple 
structure, e.g., without branching). So below we discuss some specific properties concerning 
quantum waveguides, present the effective equation of adiabatic (longitudinal) motion in 
thin films and tubes, and briefly touch upon only a few possible applications of the general 
theory to problems of quantum waves in nanotubes with spin taken into account. The results 
of this subsection represent the particular case of the general results concerning quantum 
wave propagation in thin films taking spin into account. These more general results are 
obtained together with J. Briining 

Thin crystalline films of width ~ lOnm (a few monoatomic layers), synthesized recently, 
give a more complicated example of quantum waveguide. Such a film is a waveguide for a 
quasiparticle with charge e propagating along the film, and we can affect this particle by 
means of an external electromagnetic field. In reality, a quasi-particle has spin, but we shall 
not consider spin effects for nanofilms. 

The nanofilm width do ~ Inm (10 A) is comparable with the de Broglie wavelength 
A = 27r/~ 1 nm of an electron with energy of the order of the Fermi energy ~ 1 eV. This 
circumstance leads to the following effect of “dimensional quantization” of low-dimensional 
systems: the domain of the wave function localization in the normal direction to the film has 
dimensions ~ A, and the energy corresponding to the motion in this direction is quantized. 
Therefore, the total three-dimensional problem of describing the quantum states can be 
divided into several reduced problems (on “subbands of dimensional quantization”) already 
with two-dimensional quantum effective Hamiltonians (along the him surface), which, in the 
end, allows one to obtain a sufficiently explicit description of these states by using asymptotic 
formulas. 

The him boundaries play an important role in the future constructions. A natural idea is 
to simulate the boundaries of the him by means of the Dirichlet conditions or “rigid walls” 
for the wave function. But it is more convenient to simulate by using the so-called “soft 
walls.” Boundaries of this type are related to the physical mechanism of conhnement of 
electrons near the physical him. The conhnement appears as a result of the electrostatic 
interaction between the him and the quasi-particle. One usually simulates this interaction 
by introducing a conhnement potential V\nt in the direction normal to the him. The potential 
Uint increases very fast near the imaginary boundaries of the him (the “walls”). Thus the 
wave function decays very fast outside the him and the conhnement potential Uint replaces 
the “rigid” walls simulated by the Dirichlet conditions. From the “mechanical” viewpoint, 
the conhnement potential represents the interaction with imaginary walls. The same idea 
is used in simulation of nanotubes. We shall present the corresponding formulas somewhat 
later. 

The ehective dynamics of quantum states in the approximation of the strong coupling 
method is determined by the Schrodinger equation: 

^ p2 

ih^t = '^ = n -^ Ant(r) + ^^ext(r, t), (4.22) 
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where P = —ihV — (e/c)A(r,t), e = —Cq is the charge of electron, m is the effective mass 
of quasi-particle, c is the velocity of light, and (next(i', ^), A) are the potentials of external 
electromagnetic field. We shall consider a space-uniform time-dependent magnetic field 
H = H(f). 

The characteristic value of the transverse energy e± in the tube can be found from the 
uncertainty relation: since the “transverse” momentum is ~ h/do, we have e± ~ h^/{md‘^). 
Let us introduce the characteristic “longitudinal” length Iq. Depending on the problem 
considered, Iq can be, for instance, the radius of curvature of the film or the radius of 
the solution localization area, etc. We assume that Iq 3> do- We introduce the magnetic 

length Im = \Jftcj (e|H|), the magnetic field quantum $0 = 27rhc/e, and the dimensionless 
magnetic held as the number of magnetic hux quanta passing through the characteristic area 
/odo: H' = kdo/lM^ ■ H/|H| = 27r/odo ■ H/<ho. We introduce the new variables r' = r/Zo, 
the dimensionless time t' = t/T, T = mdolo/h, the dimensionless potentials = Uint/^a, 
= '^ext/e^a, A' = edo(hc)“^A, and the dimensionless constant a' = ha/d^ and divide 
both sides of Eq. by the energy of transverse motion eo- Below, we shall omit the 

primes. Then the equation describing the motion of a quantum particle (or a quasi-particle) 
in a quasi-two-dimensional crystal takes the form 

i/iTi = = l/2(-i/iV - A)^ -1-nint(r) + next(r,f). (4.23) 


The fact that we consider the last equation in a him is determined by the boundary condi¬ 
tions. We shall assume that the him is determined by some smooth surface T. This means 
that Eq. holds and the boundary conditions are formulated in some neighborhood 

of T. It is convenient to use the special curvilinear coordinates for the description of these 
conditions, as well as for all future investigations. 

Curvilinear coordinates. By x = (a;^,a;^) we denote the (dimensionless) local coordi¬ 
nates on the surface T, then each point r in a neighborhood of T can be determined by three 
values (x^,I/), where y is the distance between the point r and its projection R{x) G T. 
Then we have 

r = R(x) -|- 2 /n(x), 

where, as above, n(a;) is a unit normal vector on T. Note that, in general, the coordinates 
x^,x^ are not orthogonal, but always (n, n) = 1, (n, djR) = 0, Z = 1,2. Thus the metric 
tensor is 


G 


ab 


7*j 

0 

0 

1 


G = detGab, a, 6 =1,2,3, 


(4.24) 


where 7 ^- = (djr, djv), 7 = det 7 jj, i,j = 1, 2, and G = 7 . Now let ns present the components 
of the vector potential A in the coordinates {x^,x‘^,y). We choose the symmetric form of 
the vector potential A = l/2[H,r]. Hence the vector potential satisfies the Lorentz gauge: 
dA°‘/dr°‘ = 0. From now on, it is convenient to use the Einstein notation and the summation 
rule. 

Soft and rigid walls. Using the curvilinear coordinates, one can consider an “empty” 
him with “rigid” walls: Uint = 0, = 0 (the Dirichlet condition) or a him with “soft” 

walls: Uint 7 ^ 0, T(a;, y) G ^ 2 ( 2 /) at each x. However, the last dehnition requires v\-at{x,y) 
to be identically dehned in the entire R^. The last condition is too strong, since T(x, y) is 
exponentially small for y ^ y and any conditions on the function T in this region ahect its 
behavior negligibly. To be dehnite, in what follows, we assume that \h(a;, j/)|aQ = 0. The 
“empty” him with “rigid” walls can be considered as the limit of soft walls described by the 
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potential rapidly increasing near the bonndary. As an example, let us consider the potential 
= {y'/D{x)Y"^, y' = y/y, m > 0 . As m —> oo, we have Vint{x,y < d{x)) 0 and 

Vint{x,y > d{x)) oo. 

Operator-valued symbol in a nanofilm. It is well known that the use of the function 
instead of the function can signihcantly simplify the corresponding calcu¬ 
lations. Substituting the function into Eq. (14. 2 dll , we obtain the following 

equation for the function 

iy% = 7^'^', n' = (4.25) 


Using the formula = G°^dad^, + G^^d^, + G-^l^da{G^I'^G°^d^{G-^l^)), we obtain 


H = iG-y^p, - - f 




(^1/4 


-G°'^AaPb —^G“^Aadfe(lnG) + -G^^AaAb + Vext{x,y,t) +Vint{x,y/y). 


We want to study solutions to Eq. that have only a few oscillations in the transverse 

direction. From the physical viewpoint, it is clear that, in general, the nontrivial behavior 
of such solutions should be determined by two-dimensional effective equations of adiabatic 
motion (HID on the surface E in a neighborhood of the physical him. Recall that our goal is to 
hnd these reduced equations corresponding to solutions with different numbers of transverse 
oscillations. As we have different scales in the transverse and longitudinal direction, it is 
natural to use the variable y' = y/y instead of y. To simplify the notation, we omit the 
prime. Then the operator Tf' in Eq. is 


'H' = '^{PiPj - ^Aipj + AiAj) + -iPl- ^Aypy + Aj) -F next(r, t) + Vir,t{x, y)- 


pj - - ^-Aydy{\n-i)- 


y 


2 7IA 


di 




1/4 


7 


1 1 


7 




7 


1/4 


(4.26) 


where Ai = (d^r. A) and Ay = (n. A). Eq. (14.251) . with the Hamiltonian determined by 
formula (lOHll . is the object of our future study. 

Using these formulas, we hnd the hrst and second terms of the expansion of the symbol 
of operator H': 


^0 + Ve^t{'R{x),t) +Vint{x,y) (4.27) 

'H'l (^x,p,y, = \yiiPiPi + (Vnext(R(2:)A),2/n)- 

+ ^Al[dy{\n-^)]y=^ , (4.28) 

where p = (pi,P 2 ), A = Aq -h yyA^i, Aq = 1/2[H,R], Ai = l/2[H,n], = A° -F yyA] -F 

0(/i2), AO = (a,R, Ao), A) = (a,R, A0 + (a,n, Ao), A, = Al +yyAl + 0{y^), A^ = (n, Ao), 
Al = (n, Ai) = 0 ^Vi=Pi- AO, and Vy = Py - A^. 

To describe the so-called slow modes, we need to compute under the assumption that 
Vi = 0, Uext = 0, and dYl/dt = 0. We obtain 


g{x) 


(xi - xa)^ 
8 


2 gV^ 


di 





(4.29) 
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This term is independent of y and py and contains only geometric characteristics of the 
embedding (the first snmmand) and the limiting manifold (the second summand). We call 
it a geometric potential. 

Effective Hamiltonians of longitudinal motion. Now we present Xo the effective 
“adiabatic” Hamiltonians The index u enumerates the Hamiltonians which, in 

our problem, is called the effective adiabatic Hamiltonian on the v-i\i subband of the size 
quantization. Substituting the function Xo = exp(ij/(n, Ao))w'^ into Eq. (Id.l2|l with Hq = 
Tfg, we obtain 

+ Vext(R(a:),t) +e'P{x), xo = exp{iy{n, Ao))w''. (4.30) 

where w'^{x,y) and s'p{x) are the respective eigenfunction and eigenvalue of the following 
problem: 

y)^ w‘'{x, y) = e^^{x)w''{x, y), w''{x, Yi{x)) = w''{x, Y 2 {x)) = 0. (4.31) 


It is well known that the spectrum of this problem is nondegenerate, thus the symbols L and 
X'" are scalar functions. For the model potential Vmt{x,y) = {y/D{x)Y'^, m > 0 considered 
in item 2, we obtain £(((x) = ((i(0)/(i(a:))^£(((0), where d{x) = D{x)^^d{0) is the dispersion 
of the state with energy s'p{x). Assuming that the width of the film is proportional to d{x), 

m 

we conclude that Zl(a;)™+i is the coefficient of homothety. As m —> oo, this coefficient tends 
to D{x). So we obtain the natural result stating that, in the model of empty him with rigid 
walls, the width of the him is equal to the distance between the walls. 

Let us present the hrst correction yLi in an expansion of the symbol of the ehective 
Hamiltonian of longitudinal motion. It is given by formula (Enni). 

Using the formula for Xo expansion of the gauge condition dpY^Aj) + 

l/ 27 *-^Ajdj(ln 7 ) + dyAy + l/2Aydy{\\ix) = 0 with respect to y = yy', it is easy to obtain the 
relations 




. 9Xo 


= iY (n, 


dAr 


OHgff dxo \ o / A \ 

^ 0 ) ^ 73 — 77 — ) =iYg ^Vidj (n, Ag), 
y 


dt /„ dt / ’ dpi dxi 

where Y = (xg, yXo)y^ and 

Yg^^V,{dyYi,AP) = -Yg^^V, (n, d,A q) = 1/2(H, A), 

S.(9"7) + \9“A«a,(\ng) + A;[S„(ln7)],.„ = 0. sfpj " A]) = S.(g«'P,), (4,32) 

From this, we hnd 


Li = -Yay’^VPPu - (E(R(a;),f),Un) - (H, A) - 


(4.33) 


where A = \Yn,V],'P = g'^^Vidj'K, and E = — Vuext — T{dAo/dt). We shall see below that 
the correction L 2 is important in the construction of the leading term of the asymptotic solu¬ 
tion to the ehective equation of adiabatic (longitudinal) motion only under the assumptions 
Vi = 0, Uext = 0, and dH/dt = 0. In this case, it coincides with the “geometric” potential 
Q{x) (I4.29p . 
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4.6 Nanophysics: wave dynamics in nanotubes 


The Pauli operator. Lengthy molecules consisting of a great many atoms situated on 
cylinder-type spatial surfaces are called nanotubes jHSl IMl IHH IHHl EH EHl EH EO] • The surface 
of such tubes can have some additional internal torsion. The nanotube diameter dg ~ 1 nm 
(10 A) is comparable with the de Broglie wavelength A = 27:/kp ~ 1 nm of an electron with 
energy of the order of the Fermi energy ep ^ 1 eV, and the nanotube characteristic length Iq 
is signihcantly larger than do- 

In the approximation of the strong coupling method, the wave functions in nanotubes 
are determined by the nonrelativistic one-particle Hamiltonian, i.e., by the Pauli operator 
with the spin-orbit interaction taken into account: 


n 


p2 

7^-K Hnt(r) + next(r,t) 


eh 

2mc 


(cr, H) -I- Tdso, 


P 


—iW — -A(r, t). 
c 


( 4 . 34 ) 


Here r G is the radius vector of a point in a neighborhood of the tube, Tdso is the 
operator of interaction of the electron spin with the electric held of the crystal EH: — 


a 


VU; 


int5 


and a is the constant of spin-orbit interaction. This Hamiltonian differs 

from that in a nanohlm only by the presence of terms describing the spin effects. 

Thus, all the notation is the same. 

In this section we consider some of results published in EH EH mi. 

Curvilinear coordinates in tubes. As in the case of thin dims, it is convenient to 
perform all arguments by using a special system of curvilinear coordinates. We assume that 
the tube axis (the curve) 7 is given by the equation r = r e M^, where R(x) is a 

smooth vector function and x G M is a natural parameter on 7 (the tube length is counted off 
from a certain point x*), |da;R(x)| = 1, dx = d/dx. If |d^R| A 0 ) if is determined the Frenet 
trihedron. The curvature k{x) = |d^R| and the torsion x(x) of the curve 7 are connected 
by the Frenet trihedron |da;R, n = d^R/|d^R|,b = [d3;R, n]} at each point x by formulas 


= — xb — kdx^i and = xn. 

By n(x) we denote the plane intersecting the tube axis at the point R(x) orthogonally 
to the axis; the section of the tube by this plane (the area in n(x)) we denote by f2(x), 
the boundary of ff(x) we denote by dVL{x). Then the tube is the union of areas ff(x), and 
its boundary is the union of dVL{x). The “physical meaning” of the boundary dVL{x) and 
of the boundary conditions will be discussed later. We introduce dimensionless coordinates 
(x, 1/1, 1/2) determined by the relations r = /oR( 3 :) + y, y = do|/ini(x) -|- (io2/2n2(x), where 
{ni(x), n2(x)} is the basis in the plane n(x). 

If we put Ill = n, 112 = b, then the coordinates thus introduced will be nonorthogonal. 
It is convenient to introduce orthogonal coordinates (see [73 EH)- First, let {ni(x), n2(x)} 
be a certain orthonormal basis in the plane n(x) smoothly depending on x (in general, this 
basis does not coincide with n, b); and let 6 *(x) be the angle between the vectors n and rii. 
Then, along with the torsion x, we can introduce an “effective torsion” Xeg = — (Sxni, 112) = 
x — dx6. Choosing the angle 6{x) (along with {ni(x), n2(x)}) so that dx6 = x, we let Xeg be 
zero. The coordinates thus constructed are orthogonal (around the tube axis, where they are 
specihed). The components of the metric tensor i,i = {x, j/i,j/2} in these coordinates 
are determined as follows: (700 = = (1 ~ gn = (722 = 1 , and 77^- = 0,i 7^ j. 

Everywhere below we shall use these coordinates. All formulas obtained below are valid in 
the case of a straight axis if we set fc(x) =0 and x(x) = 0 . If A;(x) 7^ 0 , then 771 and 772 are 
the coordinates only in the area where 1 — A:(y, n) > 0 . It follows from the considerations 
about the tube curvature given below that these coordinates are determined in the area of 
the tube axis under study. 
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Boundary conditions and geometry of nanotubes. As in the case of a nanofilm, 
the “surface” of a nanotube can be simulated by “rigid” and “soft” walls. The rigid walls 
are determined by the imaginary surface of the tube and the Dirichlet conditions on this sur¬ 
face. The soft walls are simulated by an appropriate choice of the crystal potential 
rapidly increasing while approaching the imaginary surface of the tube and creating a poten¬ 
tial well where the electron wave function is localized. Outside this well, the wave function 
is exponentially small. We shall consider tubes whose cross-section by the plane n(a;) ro¬ 
tates with respect to the basis {ni(a;), n2(x)} in which the metric tensor is diagonal, and 
simultaneously expands in the plane x = const with respect to the point R(x). We dehne 
the tube chirality as follows: we fix a cross-section Oa,* for some x* G 7 and assume that, at 
a point X ^ X*, the cross-section Oj. is obtained from Oa,* by a turn through an angle <h(x) 
(i.e., through an “angle of internal torsion” with respect to the basis {ni(x), n2(x)}) and by 
expanding by a factor D{x). 

The domain f 2 (x) can be introduced as a multiply connected domain, for example, in the 
form of a circular or an elliptic annulus. The adequacy of this representation depends on the 
form of the crystal potential in a concrete nanotube. If the domain f 2 (x) is simply connected, 
then, in the physical literature, such a nanotube is called a “quantum wire” jSlISHl. 

Operator-valued symbol. The way of introducing dimensionless variables in a nan¬ 
otube is the same as in the previous Example 4 . 5 . It is convenient to pick out the factor 
(7-1/4 wave function T, where G is the determinant of the metric tensor in the 

variables (x, 1/1,1/2), i.e., to substitute T = into the original equation. Then the 

function T' satisfies the equation = Ti'T', 7 i!' = . In what follows, we shall 

use the wave function T' and the Hamiltonian Ti'. After some transformations, the quantum 
Hamiltonian takes the standard form (Q with the operator-valued symbol 

'hL = 'H.q y.'hL'i -\- y?{Q{x') -\- 


where 

^0 = ^ + (R(x), t) + Y .]=1 {x.y), 

= -1/2(4R, H)/ ik/ 2 {n, Aq) (kiy, n)Vo - l/ 2 (y_L, H) jp-F (4.35) 

+ {|Vxext(R(2:),t) + l/2[Ao,H],y^ - l/2(cr,H) -F/i“^Q;(cr, M), 

and we introduced the notation 

Vo=p - (d^^R, Ao), Vj = -id/dyj - {nj, Aq), j = 1,2, 

i = i{y2d/dyi - yid/dy2), y± = [y, 9 ^R] = yina - 2/2ni, Aq = 1/2 R(x)], (4.36) 

M = 4R - S^V,) + 

For the symbol 7^2, we present only its “geometric” part Gix) = —k‘^/ 8 , the “remainder” 
7^2 = 7^2 — ^ is a polynomial with respect to the momentum Vo and the components of 
the magnetic held H(t) with zero constant term and with coefficients smoothly depending 
on (x, y). In what follows, we shall see that the explicit form of this “remainder” is not 
necessary for the construction of the leading term of asymptotic solutions to the effective 
equation of adiabatic (longitudinal) motion. 

The boundary conditions (rigid and soft walls) dehning the nanotube are similar to those 
in the case of a nanohlm. The corresponding change will be discussed somewhat later. 

Reduction to equations on the tube axis and the adiabatic Hamiltonian. Now 

we want to use the scheme of §3 and to hnd the symbol L of the effective equation of 
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adiabatic motion (along the tube axis). This case is characterized by the fact that the 
reduced equation contains a single spatial variable. Thus, it is natural immediately to 
separate the factor exp(i f^^(dx:R, Ao)dx/^) in the wave function. This separation takes the 
extended momentum operator V into the “short” operator p = —ifid/dx but, in the case 
of a magnetic held depending on time t, gives the correction {d^R^dAQ/dt) dx to the 
effective potential. Next, because the function T is a spinor and Tio is a scalar operator, 
the true multiplicity of degeneration of the term determining the reduced equation is equal 
to 2 r (the dehnition of r is given later). 

Taking these remarks into account, we present the solution T of Eq. (EH in the form 


T(x,|/,t,/i) = x" 


exp ( i / (c? 3,R, Ao)(ix//i ) V’ 


d 2 

y, t, /i), (4.37) 


where the symbol x''{x,p,y,t, y) = Xo{x,p,y,t, y) + pxi{x,p,y,t, p) + ... of the 
(pseudo)differential operator x’''{x,p,y,t, p) is a matrix function consisting of 2 r columns 
and 2 rows and pj is a vector function with 2 r (interacting) components satisfying 
Eq. (j3.7|l . As was mentioned above, to construct the leading term of the asymptotic solution 
to this equation, we need only to have its essential part Lq{p,x) + pL'{{p,x) + p‘^Q{x). 

Equation (ITT^ can be reduced to the equation 

+ v-mtix, y)^ = e^^{x)w'' (4.38) 


by the substitutions Xo = exp(i(y, Ao))w^ and 

^eff = y+ ^ext(R.(a;),^)+£l(a:)+ ^ (d^R{x'),^^{x\t)^ dx'. (4.39) 

Here v is just the number of the (classical) effective Hamiltonian (or the adiabatic term) 
which is also called the number of a subband of dimensional quantization. The eigenvalue 
s'f{x) is the energy of the v-t\\ transverse mode at a point x. In contrast to the case of 
nanohlms, the eigenvalues e\_ (and hence the effective Hamiltonians) can be degenerate. 
The number r, which appeared above, is precisely their multiplicity. Generally speaking, 
r can depend on x, and in this case the effect called “the intersection of terms or effective 
Hamiltonians” can occur j77l|7Hl[7SllHni. Here we assume that r is independent of x. Finally, 
we have 


\w 


1 ) 




Xo = exp(i(y, Ao) 
w^{x,y) 


w. 


E. = 


0 


1 ,... 

0 

Wi(x,2/) 




(4.40) 


wt 


{x,y) 

0 


wt 


0 

{x,y) 


where 0 is the tensor product of matrices and Eg is the unit 2x2 matrix. The matrix 
function Xi can be found from Eq. (I3.14|l . 

Now we discuss the choice of the model potential. We will consider a tube with soft walls 
and with the same elliptic cross-section IHII, which can be modeled by using the potential 


Vint{x,y) = Uint 



$(x) 

Dix) J’ 



m > 0. 


(4.41) 


Passing from the variables y = (yi,y 2 ) to the new variables y' = (yj,?/^) determined by 
the relation y = D'^y '7 = m/ (m-|-l), we obtain e’f{x) = D{x)~'^'^e’f{0). It is easy to see that 
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the dispersion d{x) with respect to the coordinates y in the state depends on x according 
to the relation d{x) = D{x)~'^d{ 0 ). It is natural to assume that d{x) is proportional to the 
linear dimensions of the tube section. Then D"< is the “soft” coefficient of extension of the 
section, and 7 is the stiffness coefficient of the walls. The dependence of the energy on x can 
be represented as 

.!(.)=.(4.42) 

As m — > cxo, the potential (ICTD disappears in the interior of the domain and tends to cx) 
outside this domain; the coefficient 7 —4 1 and d{x) D{x)~^d{ 0 ). In the limit, we obtain 
the “empty cylinder” model: nint(x*, y) = 0 for {yi/aY + (1/2/fe)^ < 1 and , y) = 00 for 

(j/i/a)^ + {y2/bY > 1 , where D{x) is the coefficient of extension (of homothety). As in the 
case of nanohlms, we introduce additional “rigid” walls in the area where the wave function 
is exponentially small. 

Taking into account the form of potential (ICTl) . we obtain the relation 

Remark. Calculating Xoi we do not £x any special form of the functions 
We assume that they form an orthonormal basis in the eigenspace of problem dCT cor¬ 
responding to the eigenvalue (term) s^j_{x) with the number v and depend smoothly on all 
its variables. Of course, such a basis is not unique, and it is convenient to make its hnal 
choice in the subsequent construction of asymptotic solutions. For example, sometimes, wf 
can be taken to be the eigenfunctions of the momentum operator I, i.e., in this case, it is 
necessary to distinguish the states inside the term according to the projections of the orbital 
momentum in these states on the tube axis. Then the momentum matrix A is diagonal. Of 
course, it is also possible to change the basis in the space of spinors; this is convenient for the 
case in which the spin affects the classical dynamics (see below the “medium-wave regime”). 
Obviously, the choice of a new basis is equivalent to the inclusion of some unitary 2r x 2r 
matrix depending on x into formula ( 1071 ) after the operator xf'. 

Effective Hamiltonians of longitudinal motion. Using formulas from § 3 , we obtain 
Lq and D). In general, the objects A, Ly, etc. introduced below also depend on the number v 
(as well as Xj) L’), and fj''). Sometimes, we omit this dependence to simplify the notation. 

The symbols Lq and Lj' are determined as follows: 

Loip,x) = HeS Er®Es, 

Ll{p,x) = ik/2{n, Ao)Er ^ E^ + Ly ^ E^ + Er ^ Lgy, Lg = -|(cr, H), 

Ly(p, x) = ({d^^)p - l/ 2 (d,R, H)) A - (Yx, H)p + (Y, Vuext + ^ + kp^^) , 

Lsy{p, x) = pL~'^a (g) (cr, da,R) ® (cr, ni) -1- (g) (cr, n2) j. 

By Ey we denote the unit r x r matrix, by A(a;) we denote the r x r momentum ma¬ 
trix with elements Ajji = , by M^{x) we denote the r x r matrix of the from 

= -i {wf, {idiVint)d2 - id2Vir,t)di)wf)^, = {wf,{d2Vint)wf)^p, = 

- {wj, {diVi-nt)wf)^p, where di = d/dyi, and by Y(a;) = Wni -F Y2n2, Y_L(a:) = Y2ni - Wna 
we denote the three-dimensional “vectors” whose components are the 2x2 “dipole” ma¬ 
trices {Yi)jji{x) = (wj,yiWj,) , i = 1,2. As above, (•, •)y denotes the integration over the 
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variables y. The symbol L 2 {p,x) is significantly more complicated, bnt we need only a part 
of it, i.e., the so-called “geometric potential” G{x) = —{k‘^{x)/8)Er (8) Eg. In the long-wave 
approximation, it is necessary to take this term into acconnt. Precisely this term generates 
bonnd states in an empty wavegnide |HHj . 

Additional boundary and initial conditions. Formulas (j 4 . 44 |l allow one to construct 
the leading term of different asymptotic solutions to the effective equation of adiabatic (lon¬ 
gitudinal) motion of i^th subband of dimensional quantization. To perform more complete 
constructions, one has to sum solutions with different numbers v. But from the physical view¬ 
point, one is interested in the reduced equations with only small numbers n, and it usually 
suffices to consider the case of several n. (In nanotubes, z/, as a rule, does not exceed 7 
This fact turns out to be very important later in the study of equations in curved waveguides 
and tubes and allows one to ignore the applicability problem for the asymptotic formulas 
obtained for large u, and the convergence problem for the corresponding series with respect 
to the number u. For this reason, it suffices to pose the additional boundary and initial condi¬ 
tions already not for the original equation, but for finite (here one-dimensional) (simplified) 
equations of the form dSIZl)- Using the physical terminology, we can say that it is of interest 
to study the longitudinal dynamics of a small set of subregions of transverse quantization. 
The corresponding additional conditions for the case of spatial waveguides will be posed 
accurately below. 

5 Asymptotic solutions to the effective equations of 
adiabatic motion 

Now we want to discuss the question about solutions to the reduced effective equations of 
adiabatic motion. The existence of the adiabatic parameter p allows one to separate the 
fast motion from the adiabatic motion (the electron motion from the nuclei motion in a 
molecule, the transverse motion from the longitudinal motion in waveguide, electron waves 
from the lattice oscillations in crystals, etc.). A very important fact is that the parameter p 
slightly depends on the energy of adiabatic motion in a certain range where it varies. The 
fact that the adiabatic approximation holds for the entire region of energies of adiabatic 
motion is well known in physical literature (see, e.g., IH2I). The adiabatic motion can be 
essentially different for different energies from this region. This fact is important if one 
is interested in the construction of asymptotic or exact solutions to the reduced equation 
of adiabatic motion describing different physical processes (and corresponding to different 
energies). This means that asymptotic and sometimes exact solutions are of different type 
and thus the process of determining the leading (or essential) part of the symbol L should be 
revised. For instance, some parts of the correction Li must be moved to the leading part of 
L. As we also mentioned, this fact can, in turn, change the definition of the characteristics of 
the reduced equation and, in particular, lead to the “semiclassical splitting” of terms in the 
degenerate case (when r ^ 1 ). In the case of nanotubes, this effect shows how the spin affects 
the determination of classical characteristics. We think that the best way to explain these 
phenomena is to consider a simple nontrivial example which, in our opinion, is the problem 
of quantum waves in nanotubes. Thus, we restrict ourselves to this example bearing in mind 
its importance. Moreover, it seems advisable to explain the main ideas and considerations 
with the example of the Schrodinger equation with the Hamiltonian 

1 

n = -Y^x- ^^y + v{x,y), ( 5 . 1 ) 
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where x G M"" and y G This is done in §5.1. If one chooses an appropriate potential 
v{x,y), this eqnation describes problems of molecnlar physics, as well as of qnantnm waveg- 
nides. The ideas from §5.1 are applied to the problem of qnantnm waves in nanotnbes in 
§5.2. 

As we mentioned above, the semiclassical analysis of the rednced effective eqnation for 
adiabatic motion is well developed and the solntions for the above-listed problems are given in 
the simplest form by nsing the Maslov canonical operator iniEsi- To obtain explicit formnlas 
is a problem which mnst be solved in concrete sitnations. There are many pnblications 
devoted to the Maslov canonical operator. Here we only note that this is actnally a certain 
algorithm whose realization, as well as the process of obtaining an answer appropriate from 
the viewpoint of the nse in applied problems (e.g., the plots of solntions, the calcnlation of 
the scattering data, the freqnency of beating, etc.), even in the one-dimensional sitnation, 
reqnires additional efforts and the nse of compnters. A detailed description of the solntions 
based on this algorithm and concrete physical resnlts are not the goal of this work. These 
resnlts present the contents of other pnblications (see, e.g., isniHni)- Here, in §5.3, we only 
very briefly describe the asymptotic solntions and the simplest physical resnlts. This remark 
also concerns all the examples considered in this paper. 

The majority of the ideas stated below in §5.1 can be generalized to other examples 
discnssed above. Nevertheless, it is necessary to emphasize that a “simple” Hamiltonian of 
the form (EH) has a very special form, and hence some effects related to this Hamiltonian do 
not occnr in the examples with other Hamiltonians. On the other hand, if the Hamiltonian 
has a different form, then other interesting effects can appear. 

5.1 General considerations. 

Internal and external parameters. In the examples nnder stndy, we implicitly assnme 
that the dimensionless coefficients (e.g., the potential next in nanofilms and nanotnbes) is 
independent of /i. Indeed, in real sitnations, it is sometimes natnral to assnme that the 
coefficients can depend on both /i and other parameters. These parameters characterize the 
kinetic energy of adiabatic motion, the strength of external helds, the strength of interactions, 
etc. In this case, the fnnctions Xg, yj', ..., Lg, Lj', ... in formnla (j3.9|) also depend on 
these parameters. Nevertheless, nnder appropriate constraints, the formnla of separation of 
variables dSH remains valid. For pnrposes of mathematical rigor, some constraints mnst 
be imposed on these parameters of the problem so as to connect them, for example, with 
the parameter /i. However, in this case, one must bear in mind that, in concrete situations, 
all these parameters are numbers, and such constraint formulas are of a very conventional 
character. For this reason, to avoid cumbersome notation, we present the explicit dependence 
on such parameters only if it is necessary. 

Semiclassical parameter h. The fact that the solutions of the equation of adiabatic 
motion can be essentially different originates from the existence of an additional parameter, 
which characterizes the excitations of the adiabatic subsystem in the allowable range. To in¬ 
troduce this parameter, we consider the hrst well-known asymptotics of Eq. (EH correspond¬ 
ing to different energies in the case of the Hamiltonian Tig = |(—+ |(~*^)^ + u)- 

Denote, for a moment, by (■, •) the inner product in the original configuration space. The 
kinetic energies of fast and adiabatic motions in this situation are Kf = (T, |(~'^^)^^) and 
Ka = (T, i(-z/x^)2T), respectively. 

As we mentioned above, according to jBj, the semiclassical solutions to Eq. (EH have the 
WKB-form T x{dS/dx,x,y)'ip{x,t, y), ip ^ exp{iS{x,t)/y)ip{x,t, y), where the function 
ip{x,t,y) depends regularly on y. In this case, the kinetic energies of fast and adiabatic 


31 


motions have the same order. Kf Ka ~ 1. This solution corresponds to the excited 
state of adiabatic motion. On the other side, Born and Oppenheimer mi constructed the 
harmonic oscillator type solution T x(a:,?/,/i)'0(a:,/i), 'ip ~ exp(—a:^//i). This solution 
corresponds to the kinetic energy ~ /x, hence the energy of adiabatic (nuclei) motion in 
this case is much smaller than the energy of fast motion. Further, in the theory of waveguides, 
we sometimes have solutions of the form T h)) where xiVy ^y lAy i^i^y ^y /^) 

depend regularly on /x. For these solutions, we obtain Ka ~ /x^. The above-listed different 
asymptotics can be classihed by the parameter h = /xy/Kf/Ka Ka/Kf ~ /x^/h^. We 
call h the semiclassical parameter. Let us emphasize that the adiabatic parameter is always 
assumed to be small and, conversely, the parameter h can be small but can also be ~ 1. 

This parameter can be explained in another way. For clarity, we consider the plain 
quantum straight waveguide and, for a while, return back to dimensional variables. We have 
the diameter do and the length Iq of the waveguide. Recall that our goal is to construct 
asymptotic solutions of the reduced equation which describe the motion along the tube axis 
in a sufficiently wide range of longitudinal energies and the transverse wavelength A_l ~ do- 
To the longitudinal energy, there corresponds the characteristic de Broglie wavelength Ay = 
h/p||, where p\\ is the dimensional momentum of longitudinal motion. Now the “semiclassical” 
parameter is h = A||/Zo- In other words, the parameter h determines the “smoothness” of 
the function ip {h~^ is the number of oscillations at the distance ~ Iq) and agrees with the 
estimation of its derivatives: ("0, |^) ~ h~^. 

We again consider the example of an empty waveguide. Then the energy of the longitu¬ 
dinal motion on the ix-th subband of transversal quantization has the form: p‘j/2m + Veff{x), 
where Ves{x) = Uext(a^) + r‘^7i‘^s±, e± = h^/{2mdl). Denote by p± the transverse momentum. 
Taking into account the relation between the de Broglie wavelength and the corresponding 
momentum, we obtain ^ Thus, the kinetic energies of longitudinal 

(adiabatic) and transverse (fast) motions satisfy the relations Ka/Kf ~ pP^/h? and Kf ~ e±. 
Now we return to dimensionless variables. Then the dimensionless longitudinal kinetic en¬ 
ergy is Ka = p|/2 ~ p?/h?. It is clear that if a particle moves along the waveguide, then 
the kinetic energy can vary under the action of the force / = —dves/dx. For this force / 
not to accelerate the particle so that its kinetic energy be of order different from p‘^/h‘^, 
it is necessary that its work does not exceed, in the order of magnitude, the parameter of 
the characteristic kinetic energy corresponding to the initial momentum. In dimensionless 
variables, to the distance ~ /q there corresponds an interval ~ 1. Hence the work of the force 
/ is of the order of the derivative dve^/dx. This implies that the effective potential must 
have the form Ues = + ^'^esi.^)y where = const and u/ff(a;) can, in general, regularly 

depend on the parameters p and h. Moreover, the work can be even equal to zero, since the 
characteristic longitudinal momentum is determined not only by the variable part u/g(x) of 
the effective potential, but also by the “input” momentum of the wave packet under study 
(i.e., by the gradient of the phase of the wave function at the initial time instant in the 
Cauchy problem or by the momentum of the incident wave in the scattering problem). In 
the last case, the asymptotics of the wave function can be obtained by using the well-known 
Born approximation. 

Remark. We point out that Ugfr is determined by both the external field and the field 
of the crystal. Therefore in the case of a quantum waveguide, the above constraints lead, 
in particular, to the assumption that the geometric parameters of the waveguide, i.e., the 
curvature (and torsion) of its axis, the width, etc., vary sufficiently “slowly.” 

Characteristic time scale and the reduced equation consistent with this scale. 
The question concerning the time scaling is nontrivial and, generally speaking, can be solved 
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separately in each concrete problem. It is natural to understand what characterisitc time is 
required for a quantum particle to walk through a certain characteristic distance. For prob¬ 
lems in nanotubes, the characteristic distance is the total tube length (e.g., in the scattering 
problem or in the problem of the wavetrain propagation) or the size of the “localization area” 
of the wave function in the problem of bound states. For a while, we assumee the character¬ 
istic distance to be of the order of /q in dimensional variables or to be ~ 1 in dimensionless 
variables. One has to replace the time scaling by the energy scaling in the case of stationary 
problems (for instance, in problems of electron-phonon interaction or in molecular physics). 

Remark. To introduce the characteristic time scale in the general case, we can use 
the following ideas. It is clear that the characteristic time scale is t ~ a/ (v), where a is the 
characteristic distance for adiabatic motion and (v) is the mean velocity. Generally speaking, 
a as well as (v) depend on the “longitudinal” kinetic energy. In quantum mechanics, we have 
(v) = d{x)/dt. Using Ea. ()2.2|l . we obtain d{x)/dt ~ ifi ^([p^/2,x]) = (p) = p/h (cf. Remark 
3 in the 113. 3 j) . Thus we have t ~ {h/fi)a. For the scattering problem, the wavetrain 
propagation problem and some other problems, we can set a ~ 1. For lower bound states 
and trapped modes we have to set a 

The dimensionless time t used in the general scheme for Eq. (HI was actually chosen 
for the case p± ~ p\\, i.e., for the case /r/h ~ 1. If this relation does not hold, then the 
time of passage of a particle through the waveguide, which is naturally understood as the 
characteristic time of the problem, must be multiplied by the factor {h/p). Therefore, instead 
of t, it is convenient to introduce a new dimensionless time t' by the relation t = {h/p)f. 

In the case of nanotubes this redefining of the time scale becomes consistent with the 
preceding physical argument because of the following transformations in Eq. (ITTIl . The 
term results only in a displacement (renormalization) of the energy in the stationary 
problem generated by the reduced equation; or the factor exp{—iv^gt'/p) appears in the wave 
function of the nonstationary Schrodinger equation dUD- Taking this into account, we 
represent the solution of this equation in the form = exp{—iv^Qt'/p)'ijj"^, where iIj"^{x, t) is 
a new unknown function. Since we assume that p ~ p/h, it is natural to divide the equation 
by the parameter p^jh?. In the left-hand side this gives the derivatives i{h‘^/p)^^, which, 
after the above change of time, take the form ih^^. It is important to point out that this 
transformation concerns only the time variable: the variables x and y are not transformed. 
As a result, instead of Eq. (EH) and Eq. (ICTll taking the corrections L^, ^^c., into 

account, we obtain the equation (the primes are omitted): 
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Accuracy of asymptotic expansions. The number of terms in the expansion of the 
intertwining operator y and the operator L can be arbitrarily large. However, as we already 
mentioned, to calculate terms of these series explicitly, even terms with small numbers, is, 
as a rule, a very complicated problem. Therefore, it is natural to consider only the terms for 
which one can correctly estimate the leading term of the asymptotics of the wave function or 
of the energy value. It is reasonable that the notion of the “leading” term of an asymptotics 
can be determined not only by the adiabatic parameter p, but also by the “semiclassical 
parameter” h, which is related to the form of the coefficients and the solution of the effective 
equation of adiabatic motion. We shall turn back to the question about numbers of terms in 
the intertwining operator y and the operator L later. Now we recall well known estimates 
which allow us to estimate these numbers. 
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Taking in account this fact, let us discuss the problem of choosing the number of terms 
in the expansion of the symbols of the operator L and the intertwining operator y. Again 
we restrict our consideration to the case of nanotubes, although the main ideas can be 
generalized for the majority of the adiabatic problems listed above, including the non-self- 
adjoint problems like water waves in a picnocline (problems of such type usually appear in 
hydrodynamics). 

Since the problem contains two parameters and h, we shall calculate as many terms 
as we need to construct the leading term of the asymptotics with respect to max(h, /i) if 
h -C 1 and with respect to /r if h = 1. (Recall that the parameter fi is always assumed to 
be small, and the parameter h can be either small or of order 0(1).) To find the minimal 
reasonable number of terms in asymptotic expansions, it is natural to use the well-known 
estimate for the solution of the Cauchy problem for a nonhomogeneous Schrodinger type 
equation: ie^ = A(j) -1- /, (j)\t=o = 0. Here A{t) is an essentially self-adjoint (for each t) 
operator in L 2 , and e > 0. Let 0 be a solution of this problem; then the following inequality 
holds for any t from the fixed interval [0,T]: ||0 ||l 2 < ^max^gjo^T] II/IU2- 

We assume that V'ex is an exact solution of the original equation (EZD and Tas is its 
asymptotic solution of the form (O and these solutions coincide at zero instant; moreover, 
0as satisfies the original equation with discrepancy /as- For the operator A, we choose the 
original quantum Hamiltonian and set e = fx and T = ^Tq, where Tq is independent of h 
and /i. Then we obtain the estimate 

h 

||0as-06x11 < ^Tomaxll/asllij (5.3) 

for the difference 0 = 0as ~ 0ex. This implies that the minimal reasonable number of terms 
in the expansion of the symbols of the operators x L in formulas (Eni) and (EHl) must 
at least satisfy the condition ^||/as 11x2 -C 1 for /i 1. Of course, it should be remembered 
that the norm of the discrepancy /as depends on /x and h. As heuristic arguments, it is also 
useful to apply the estimate o to the reduced equation (EZD. 

Classification of quantum states for longitudinal motion We return to the passage 
from Eq. (13.7|) to Eq. (15.2|) . For h <C 1, to construct a wave function, it is natural to use 
the semiclassical approximation. Outside a neighborhood of the focal points (the turning 
points), the typical asymptotics of a wave function with characteristic wavelength Ay ~ h is 
given by the WKB-solution 

0(a;,f) = A(a;,f, h) exp , A{x,t, h) = A{x,t,0) + 0{h), (5.4) 

where S{x,t) is the phase and A{x,t,h) is, in general, the vector amplitude. As it is 
known [H] , in the first approximation, after the substitution of this function into the original 
equation, the operator -ih-^ is, in the leading part, replaced by and thus the order of 

. 2 ^ 

the terms in the operators (x, f (—fh^)) in dSH is determined by the order of the 

functions Lj{x, ^|f). This fact leads to the well-known conclusion that the phase S{x,t) is 
determined by the classical Hamiltonian system whose Hamiltonian is the leading part of 
the symbol expansion with respect to the parameter h. Bearing this in mind, it is natural 
to define the operator = —ihd/dx, rather than the operator p = —ipd/dx, to be the 
momentum operator. Clearly, for p = h, the classical Hamiltonian is the effective Hamilto¬ 
nian dsm, but if the adiabatic and semiclassical parameters /x and h are of different orders 
and h -C 1, then to construct semiclassical asymptotic it is required to write the expansion 
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with respect to parameter h, assigning /i = As we shall see later, in some cases, addi¬ 

tional terms from L'( will be inclnded in the classical Hamiltonian (snbject to the expansion 
with respect to h). 

Now let us discuss how many terms in the expansion of the operator in the right-hand side 
of (j5.2|l we must have to hnd the leading part of its asymptotic solution. By setting e = h and 
applying the estimate (j5.3|l to Eq. (j5.2p . we see that, at least intuitively, it suffices to calculate 
the effective Hamiltonian (p^)^/2 -|- and the first correction ^p^)- This conclusion 
is consistent with the well-known fact from the theory of semiclassical asymptotics: terms 
of the order of (and even of 5 > 0) do not affect the phase S{x,t) and the leading 
part of the amplitude A{x,t,0). This conclusion has a general character and holds always 
if it is assumed that /x < h -C 1, i.e., for the case in which the semiclassical approximation 
can be used. At the same time, as we just noted, concrete formulas can essentially differ in 
the following situations: a) if p and h have the same order {p ~ h) and b) if /x h. If 
the parameter h ~ 1, then the semiclassical approximation cannot be used, but Eq. (15. 2 j) 
remains to hold and can even be simplified, although in this case a part of L 2 must be 
taken into account. The existence of these differences results in the following classification 
of asymptotic solutions depending on the relation between the parameters /x and h (or, 
which is equivalent, depending on the relation between the longitudinal and transverse 
Kj kinetic energies in the waveguide). 

a) For h = /X, we have the standard “semiclassical” situation [Hj or the “short-wave” 
regime in which the “longitudinal” energy is of the same order as the energy of transverse 
motion and d ~ Ay Iq in dimensional variables. Then the effective adiabatic and semiclas¬ 
sical Hamiltonians, as well as the corrections to them, coincide, and to find the leading term 
of the asymptotics of the wave function, the complete description of the effective Hamiltonian 
and the first correction is required. 

b) In this case, which is naturally called the “medium-wave” regime, /x -C h 1, the 
“longitudinal” energy of the mode is significantly less than that of the “transverse” mode 
and d < A|| < v^/od in dimensional variables. Then, expanding the correction L\ with respect 
to the parameter p/h, for the symbol of the operator, we obtain 

u = (htA + 0)) + 0)p'‘ + A ■ 

This implies that the nondifferential part ^L\{x,t]) of the first adiabatic correction in the 

expansion of the operator L\ can be transferred into the semiclassical effective Hamiltonian. 
This is clearly seen in the case h = Then the semiclassical effective Hamiltonian 

becomes equal to -f L\{x, 0). Moreover, for p, the term ^L\{x, 0) can play 

the determining role. Then an argument similar to that in item 2.6 shows that this term 
“accelerates” the particle in the longitudinal direction so that the characteristic longitudinal 
momentum in dimensionless variables takes the value ^JJi. In other words, in this case, for 
the parameter h we must take the parameter ^//x, and we return to the situation considered 
above, but with multiplied by a small value. Clearly, if L\{x,0) = 0, then the above 
argument is meaningless. But, as we shall see later, such a term appears in nanotubes both 
due to their geometry and due to the external electromagnetic field. In this case, there 
arise some additional parameters, e.g., field amplitudes, and these parameters can effectively 
decrease the value of Lj' (x, 0) and thus compensate the increase caused by the parameter 
h?I /X. We also note that the terms in the operator L\ containing the second- and higher- 
order derivatives (corresponding to higher powers in the expansion of Lj' with respect to 
the variable p^) can be omitted in the calculations of the leading term of the semiclassical 
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asymptotic solution, although, of course, these terms do not decrease the accuracy. 

Remark. If the main part of the Hamiltonian is quadratic with respect to the momenta, 
the Hamiltonian preserves its form. In the other examples considered above (e.g., for equa¬ 
tions with rapidly oscillating coefficients, for waves in picnocline), the leading part of the 
adiabatic Hamiltonian Lq depends on the momentum p. Therefore, this expansion changes 
the structure of the leading part. For example, in the problem with rapidly oscillating co¬ 
efficients (electron waves in crystals), Lq is replaced by its expansion with respect to p in 
a more complex way (non quadratically). Moreover, this expansion, as a rule, begins with 
terms quadratic with respect to p. The coefficient before p^ is inversely proportional to the 
effective mass (see jl]). 

c) If the parameter h ~ 1, then the semi classical approximation cannot be used, and 
the wave functions oscillate, if at all, rather slowly and in the dimensional variables Ay ~ /q. 
According to the above, this situation is possible only if Li{x,0) = 0. But the adiabatic 
approximation works, and from (15.21) one can easily derive the equation for the leading term 
of the asymptotics of the (smooth) wave function on the waveguide axis. For this, 

it suffices to set h = 1 in (EH and then to let /r —>• 0. As the result of this passage to the 
limit, which is naturally called the long-wave approximation, we obtain the equation 



■f 

2 dx‘^ 




.dLl 

dp 


(a:,0) 


A 

dx 


+ 0 ) 


tf. 


Remark. As already noted for the Helmholtz operator in plane one-mode waveguides, 
such an equation was first obtained in isa. where, in particular, it was proved that one 
can organize a single bound state in the waveguide by choosing an appropriate curvature of 
the waveguide. An equation similar to and several consequences of it were obtained 

in jSH E21 EBl EH EH EH EH EH] • We also note that the equations of this type are close to 
the equations obtained as the result of averaging or homogenization on the theory of wave 
process in the media with rapidly varying characteristics EOIIHSIEZIEHI. 

d) Finally, we can consider the case in which p ^ h or, in dimensional variables. Ay -C 
A± Iq- This case is naturally called the “ultrashort-wave” regime. For Eq. (Q to be 
meaningful, it is necessary to impose additional constraints on the behavior of the functions 
Lj in the variable p. We consider only the case for which LJ are polynomials of a degree not 
exceeding 2 with respect to p. Then it is easy to show that the semiclassical approximation 
for any p ^ h can be applied to Eq. (15.2|) . However, this is not sufficient for reconstructing 
the asymptotic solution of the original Schrodinger equation in the waveguide from the 
function ifi{x,t) by formula (jd.bj) . For example, if the first correction Xi in fhe expansion of 
the symbol of the intertwining operator in formula (EH) depends linearly on p, p > hf, 
and a rapidly oscillating function of the form (E3D is taken to be if, then the function pxi4> 
turns out to be not small. In this sense, the expansion of the operator x in powers of p is 
not an asymptotic expansion. For this reason, as we shall see below, the ultrashort-wave 
approximation for curved nanotubes can be used for the case in which p h. We 

can also note that, in the ultrashort-wave case, the actual effective potential is small so 
that it can be transferred into or even omitted at all. Then the semiclassical effective 
Hamiltonian coincides with the Hamiltonian of a free particle, and hence, in this case, the 
semiclassics is simply the Born approximation. 

Remark. It should be remembered that, in concrete calculations, the above classifi¬ 
cation (in the parameters p and h) must be made more precise, which is related to the 
values of both the external fields and the crystal field, as well as the relations between 
them. Of course, in this case, the corrections can also be included into the leading part 
of the symbol (in the classical Hamiltonian), which, however, can unnecessarily complicate 
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the procedure for constructing the asymptotic solutions. On the other hand, as already 
noted, in a real situation, each parameter is a concrete small number, and hence the fur¬ 
ther detailing of how the effective Hamiltonian and corrections to it depend on the relation 
between the parameters /x and h has an academic, rather than practical, character. Tak¬ 
ing this consideration into account, it is convenient, from the mathematical viewpoint, to 
hx the relations between the parameters /i and h, assuming for the respective regimes that 
a) h = /X, h)h = y/Jl, c) h = 1, d) h = /x^/^ and including the “remaining” parts of 
the relations between /x and h into the coefficients of the equation (such as the strengths of 
electric and magnetic helds, curvature, etc.). 

The suggested classihcation can be used in the general situation. But in some cases 
some regimes (like ultrashort wave region or long wave region) does not exist. Anyway the 
question of existence of solutions to a certain hxed h should be solved individually. 

The number of terms in the intertwining operator. We assume that we use the 
dimensionless time t consistent with the semiclassical parameter h (see subsection Charac¬ 
teristic time scale and the reduced equation consistent with this scale). Suppose 
that we have preserved N terms in the expansion of the symbols of the operators x'" and L'^ 
determined by the coefficients x’j and L^. Suppose also that we have constructed a function 

satisfying the reduced equation Q with accuracy up to a discrepancy /. It follows 
from the formulas of §3 that, for appropriately dehned x'j and L^, the substitution of the 
function T = x'^'ijj^ into the original equation gives 

+ rf- 

ot /x^ 

Here F is, in general, a pseudo differential operator that does not change (as well as the 
operator x'') the order of the functions with respect to the parameter h if the functions 
oscillate with the characteristic wavelength not less than h. Applying the estimate (15.3|) to 
this equation, we readily obtain the following conclusion: the function T differs from the 
exact solution by a value of the order 0{h^), 5 > 0, for h -C 1 or of the order 0(/x) for h = 1 
if A^ > 1 and the discrepancy / is equal to for h -C 1 or to 0(/x) for h = 1. Thus, 

the minimal reasonable number of terms in the expansion of the operator x^ in constructing 
the semiclassical asymptotics is equal to 2 (i.e., we must consider the zeroth- and hrst-order 
terms). But if we are interested in the long-wave approximation (i.e., in the case h = 1), 
then, obviously, in the expansion of x'^ we must consider three terms of the expansion (i.e., 
N = 2). In this case, it suffices to solve the reduced equation (15. 2 j) up to 0(/i). 

The problem of calculating the symbols of the operators and is close to the problems 
in perturbation theory for operators with discrete spectrum (in particular, matrices), and 
the function L’-' is similar to an eigenvalue, while is similar to an eigenfunction. The 
terms of the expansion of the symbols x^ and are calculated successively, but the explicit 
calculation of Lp precedes the calculation of x'j and is based on the fact that x’j exists. 
On the other hand, the leading term of the asymptotics is already determined by Xo (and, 
naturally, by the function Thus, in the construction of the semiclassical asymptotics 
in the minimal reasonable approximation, explicit formulas are required only for Lq, and 
Xq, while in the construction of the long-wave asymptotics, explicit formulas are required for 
Lg, and Xo, Xi■ Moreover, as was already discussed in item 2.9, to obtain the medium- 

wave and long-wave approximations (h S> /x), for L^, it suffices to calculate this function and 
its hrst-order derivative for p = 0, while for L^, it suffices to calculate this function for p = 0. 
This fact is very important for deriving explicit formulas. It should be noted that, although 
the terms px\'ip'^ for h -C 1 and (/xXi + h^X2)l^^ h ~ 1 are only corrections to the leading 
term Xol/’^ of the asymptotic expansion, it can be described correctly only if the existence 
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of such corrections is guaranteed, while the discrepancy obtained by the direct substitution 
of the function Xo''P'^ iiifo the original equation is not sufficient to prove that the difference 
between Xo''P'^ cind the exact solution is small. Here we have the following distinction from 
the standard semiclassical approximation for the scalar Schrodinger equation; in the latter 
case, the substitution of the leading term of the asymptotics into the original equation readily 
gives the desired small discrepancy 0{h?). 

Essential part of the quantum Hamiltonian for longitudinal motion. Thus, 
it follows from the analysis performed in the preceding items that, in all cases a)-d), the 
leading term of the asymptotic expansion of (as well as of T) is completely determined 
(sometimes, with excessive accuracy) by the quantum Hamiltonian 


t = - {-ih^ 
2 \ dx 


1 

+ ■^60“^ 

/i 




This Hamiltonian is naturally called the essential part of the Hamiltonian of longitudinal 
motion on the n-th subband of size quantization. 

Semiclassical splitting of adiabatic terms. For h 3> /x, the correction Li(a;, 0) 
can change the effective Hamiltonians. This change can be essentially important for the 
case in which the term is degenerate or (and) the original problem is a vector problem 
(i.e., the original quantum Hamiltonian is a matrix Hamiltonian). Let us consider this 
situation in more detail. Because the original problem is self-adjoint, the matrix Li{x, 0) is 
Hermitian. We assume that all its eigenvalues Ai(a;),..., \r{x) are distinct and, along with 
the eigenvectors, smoothly depend on x. Moreover, to avoid any renormalization of energy, 
for simplicity, we assume that \j {x) varies along the waveguide and \j = 0 at the beginning 
of the waveguide. Since any adiabatic term is degenerate, there is a great ambiguity in the 
choice of vector functions xfiJ = any variation in their choice naturally leads to 

a change of the matrix Lj'. We choose the vector functions yj, j = 1,..., r, so that the matrix 
L'^{x, 0) be diagonal. If 1 3> h 3> /i (we have the “medium-wave” approximation), then the 
expressions ^\j{x) must be added to the semiclassical Hamiltonian. Following an argument 
similar to that in item 2.6, we conclude that the parameter h must satisfy the inequality 
^/JI. > h 3> /r. We see that the semiclassical Hamiltonians (terms) + ^1^ +^\j become 

different for different j and thus the “semiclassical” separation of “adiabatic” terms occurs. 
The value of this splitting depends on the number ^ and attains its maximum at h y/I. 

In the hrst approximation with respect to the parameter the states corresponding to 
different semiclassical Hamiltonians do not interact. Thus, inside an adiabatic term, we can 
asymptotically diagonalize the system for the vector function ip''. Although this asymptotic 
fact turns out to be valid for any values of h from the above interval, it is clear that, for 
values of h close to /i, the subsequent terms of the expansion can be large so that the above 
representation of such an asymptotic “diagonalization” becomes meaningless. Therefore, the 
situation with h close to /i must be considered as a situation with a degenerate term, and it 
is more natural to leave the term in the original correction which allows one more 
adequately to take into account the interaction between the states inside an adiabatic term. 
This reasoning concerns the case of the matrix correction L^. In the scalar situation, the 
term ^Lj'(a;, 0) can always be included into the effective Hamiltonian, but, in general, this 
does not hold in the vector case. We have no opportunity to discuss the problems touched 
upon in this section in more detail. We only point out that all phenomena mentioned in this 
and in the preceding items are completely determined by the essential part of the matrix 
analog of the Hamiltonian of longitudinal motion on the z/-th subband of size quantization. 
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5.2 Essential parts of the effective Hamiltonians for longitudinal 
motions in nanotubes for different regimes. 


The goal of this subsection is to realize the ideas from the previous subsection concerning 
the problem of quantum waves in nanotubes with spin taken into account. 

As already noted, the “rapidity” of the longitudinal mode is determined by the parame¬ 
ter h. For h 1, the corresponding states (in particular, the levels of longitudinal motion) 
are described by semiclassical asymptotics. Since the strong electric or magnetic held in¬ 
creases the longitudinal energy of a particle while it passes through the tube, the parameter h 
cannot be chosen arbitrarily and must be consistent with the strengths of the external helds 
contained in the problem. For mathematical rigor, we must assume that there is a func¬ 
tional dependence between the parameters /i, h and the strengths of the helds. From the 
physical viewpoint, it is natural to speak about the value of diherent terms in (j4.44|l in 
dependence on the “rapidity” of the longitudinal mode and the strengths of external helds. 
For a solution with wavelength ~ h to exist, it is necessary that the work of the ehective 
held Ees{x) = —dvext/dx — de±{x)/dx— {d^H, dAo/dt) along the tube do not exceed, in the 
order of magnitude, the characteristic kinetic energy of the longitudinal motion: e\\ 

This implies the following constraints on the external helds and the “oscillation” in the di¬ 
mensions of the section. We assume that the potential of the external electric held is equal 
to zero at one of the tube ends; thus, Uext(R-(a;), t) will be equal to the work of the held 
along the tube. This implies that the external potential in dimensionless variables has the 
order of jj?/h?. We introduce the functions vl^^{x,t) and \''{x) determined by the relations 
Uext(R'(a^), ^±ix) = £([(0)-|-(/i^/h^)A^(a:) and we assume that vl^^{x,t), 

X'^{x) take values that do not exceed unity with respect to the parameter. 

In what follows, for simplicity, we assume that the function in fact, smoothly depends 
on t', and we write vl^^{x,t'). The constraint on the ehective held implies a constraint on 
the character of the time-dependence of the magnetic held. We assume that H = Hq + 

We replace the “adiabatic” momentum operator^ = —ifid/dx by the “dynamic” operator 

= —ihd/dx, obtain p = {fi/h)p^, and then divide the adiabatic Hamiltonian by p?/h? 
to “compensate” the redehned momentum. This leads to the redehned classical momentum 
p = {p/h)p^ and, by (I4.44|l . to the following formula for the symbol C'^{p^,x,t', p, h) of the 
essential part of the quantum effective Hamiltonian of the longitudinal motion: 


C = 


h\2 


(P ) 


+ -1- (j){x, t') + y{x) 


Er ® E,+ 


p 


1 1 

-Er 0 -(CT, H) - -(^xR, H)A ^Es + Lsy 


+ h 


(4<F)A-(Yx,H) 


Esp'^E 


+r (y, vti. + ^ + H/rn'j - 


h‘^k{x) 


Er 0 Eg, 


(5.5) 


where Ag = (1/2)[R, Hi(F)], 0(x, F) = (1/2) {dxll{x'),[ll{x'),dlii{t')/dt']) dx\ Omit¬ 

ting the primes, we obtain the desired reduced equation on the subregion of dimensional 
quantization in the form: 


ih^ = C’^ij’^, = C{-ih—,xH,p,h). (5.6) 

Along with formulas (1071) and (lonii . this equation, for diherent relations between p and h, 
determines the leading term of the asymptotic solutions of the original equation (El). Some- 
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times, several terms in (ESI) can be omitted, and several terms from the “adiabatic” correc¬ 
tion (sometimes, matrix terms) can be transferred into the leading part of the “semiclassical” 
effective Hamiltonian determining, in particular, the dynamics of the classical motion of a 
particle in a thin tube or, if the terminology introduced in HSIEHI is used, determining the 
nonstandard characteristics of the original equation (ESI- We describe this “transfer” and 
the corresponding classical systems in the next subsection. 

Equations of classical mechanics in nanotubes with spin and the term multi¬ 
plicity taken into account. Below we perform different expansions of the function (15.bj) . 
For h -C 1, we set = £q -|- hC\ where Cq denotes the terms larger than h and 

hCi denotes the terms ~ h. According to |15l 138 j . these terms allow one to reconstruct the 
leading term of the semiclassical asymptotics completely for h -C 1, which outside the focal 
points has the form of WKB-solutions ~ exp{iS^{x,t)/h)A^{x,t). Their phases S^{x,t) 
can be determined by integrating the one-dimensional Hamiltonian system 

/ = dHipg/dx, X = -dH^s/dp’^ (5.7) 


with the classical Hamiltonian i7gg(p^, x), which is an eigenvalue of the matrix symbol C^. If 
the matrix £g has distinct eigenvalues, then the semiclassical splitting of the adiabatic term 
occurs, i.e., several distinct classical Hamiltonians can correspond to the same adiabatic 
term Ffeff- The vector part of the asymptotics can be found from the linear “polarization” 
equation which contains C\. 

We shall consider the following situations corresponding to different relations between 
the parameters p and h. 

a) Short-wave regime: h = p. Then 


C’P = H^sEr^E,, = 


dAj 


eff = = vLti'^ix),t) + (p{x, t) + X'^{x), 


£(( = ( Y, Vuext + -^ + [(d,,<F)A - (Y^, H) 

-Er ® ^(CT, H) - H)A <^Es + L^y 


EsP^+ 


(5.8) 


The Hamiltonian system in this case is equivalent to the Newton system x = —dv^ff/dx. By 
the estimates in items 2.8, 2.10, the leading term of the semiclassical asymptotics, which is 
determined by these classical equations, must give a good approximation in problems with 
external helds |H| < 1 T, E^xt = iVuextl ^ 10“^V/nm. 
b) Medium-wave regime: h = In this case, we have 


py _ 
^0 — 


+ vlxt(R{x),t) + Hx, t) + 


w = 


-Er (8) ^(cr, H) - ^(d 2 ;R, H)A ® + Lsy 


Er®Es + W, (5.9) 

, ^=l(d,^)A-(Y^,H})0Ey. 


Thus, here the symbol £q is a matrix symbol. Here the classical effective Hamiltonians are 
the eigenvalues of the matrix £q. Obviously, they can be represented as the sum of the 
function (p^)^/2 v^^^(R(x),t) + <P(x, t) -|- A'^(x) and the eigenvalues of the matrix W which, 

in general, depend on p^. If they are distinct for all (p^, x), then the semiclassical splitting of 
the adiabatic term occurs. Obviously, these classical Hamiltonians depend on the spin terms; 
hence the spin of a particle affects its phase trajectory as follows: the energy of longitudinal 
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motion in the tube, in contrast to the total energy and the energy of transverse motion, is 
comparable with the spin energy. 

The eigenvalues of the matrix W depend on the form of the transverse section of the 
tube and on the external fields. As the simplest example, we consider the case of the 
cylindrically symmetric crystal potential (yi = gcos(f),y 2 = ^sin0) Ujnt = 'yint(^), and w’’' = 
exp{±ii>(j))u‘'{g). In this case, Y = 0, the eigenvalues of the matrix A are ±z/, = 0, 

is the diagonal matrix with elements Tkz/, k = 27r f^(dvint(g)/dg)(u‘'(g))^ dg, and 
the adiabatic term splits into four semiclassical ones determined by the classical effective 
Hamiltonians (^feff)u ^ + (^efr)n) 


(^'eff)n = t) + X’^ix) T -(^xR, H)z/ + CTti 


HzLufi 


, (5.10) 


where = ±1. In this case, we must assume that |(1/2)H ± /i“^aKz/(9a,R| ^ 0. Otherwise, 
the effect of intersection of terms or of a change in the multiplicity of the characteristics 
occurs, where the standard semiclassical approximation does not work (see, e.g., j3]). In this 
example, the Hamiltonian system is also equivalent to the Newton system with the effective 
potential (On- The semiclassical approximation must work for external helds |H| < IT, 
Eext = |Vuext| 10“^V/nm. We also note that Eq. (I5.0|l can be treated as an equation 
with operator-valued symbol (its operator symbol is a matrix), and we can again apply the 
“operator” reduction to this equation. 

c) Long-wave regime: h = 1. In contrast to the above regimes, for the “long-wave” 
approximation to exist, it is necessary to impose additional constraints on the value of the 
magnetic field. Formally, this follows from the existence of a term ~ /y in the Hamil¬ 
tonian dn^. The physical cause consists in the following. Even without the spin effects 
taken into account, the longitudinal and transverse modes are related due to the interaction 
between the magnetic field and the transverse orbital momentum, which is described by the 
term —(2/i)“^(9j;R, H)A ® Eg. In this case, the transverse energy signihcantly exceeds the 
longitudinal energy; hence the dynamics of longitudinal motion becomes very sensitive to 
small variations in the transverse energy. In turn, the transverse energy changes because of 
the variation in the magnetic held hux through the tube section due to a change in the angle 
between the plane of the tube transverse section and the vector of the magnetic held H. For 
the magnetic helds ~ 1T chosen above, this interaction accelerates the particle to energies 
that are not consistent with the parameter h. To avoid the acceleration, we must take into 
account that = £(((0) -|- /i^A^(a;), and Hq = 0, H = /iHi. This means that we con¬ 

sider magnetic helds < 10“^ T. We also assume that the dimensionless constant of spin-orbit 
interaction is a ~ /i^. This implies that the leading part of the Hamiltonian contains only 
the terms L2(0, a;)|H=o, •!)ext=o, a=o- Subsequent calculations of this term (see item A3 in the 
Appendix) lead to the last term in (15.5j) . As a result, we obtain the limit reduced equation 
of the form 


1 


2 dx‘^ 


+ t)^ Er^ Eg + (c?x*h)A 0 Eg + R" 


'ip'' = 0 


W = ( X’"{x) - y'j Er®Eg- ^(93 ,R, Hi)A ®Eg-Er® ^(cr. Hi) -F y ^Lgy. (5.11) 


Here the semiclassical approximation cannot be used and any information about the solutions 
of the reduced equation (more precisely, about the systems of equations) can be obtained 
either by general qualitative methods or by numerical integration ISH E21 ESI EH IS3 ESI 
EZlESlEi. As already noted for the Helmholtz operator in plane one-mode waveguides. 
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such an equation was first obtained in isa. where, in particular, it was proved that one 
can organize a single bound state in the waveguide by choosing an appropriate curvature of 
the waveguide. An equation similar to and several consequences of it were obtained 

in IHH] , Eq. ()5.11|1 works in the case where the correction is small. This imposes 
several constraints on the value of the constant-in-time component of the magnetic held 
Ho ~ /i and on the constant of the spin-orbit interaction a ~ /i^. 

d) Ultrashort-wave regime: h fi. In the case of ultrashort-wave approxima¬ 

tion, the external magnetic and electric helds can be stronger than those considered above. 
Namely, Uext ~ ~ /r/h, |H| < 10 T. Introducing the notation Uext = 

H(t) = (/i/h)Hi(t) and omitting the primes, we obtain 

^ Er®Es + (h + /i (^-(Yx, H)/ + (^Y, Vuext + ^ + 

El = (p^4<h) - 1/2(4R, Hi)) A - 1/2(0-, Hi) + (g) (o-, ni)p'* + ® (o-, na)^'^), 

(5.12) 

where ifgg = (p ^)^/2 -|- Ve^t{x) -1- </>. Although £q is a matrix symbol in this case, its leading 
part is Ej.®Es. Nevertheless, the matrix correction fi > h can lead to the splitting of 
the adiabatic term and, in particular, to additional terms in the phase of the wave function 
~ /r/h, /i^/h, ..., which are, in general, different for each of the states contained inside the 
adiabatic term. 



5.3 Examples of asymptotic solutions for some problems in nan¬ 
otubes 

As we mentioned above, here we have no opportunity to describe in detail asymptotic solu¬ 
tions of the effective equation of adiabatic motion in nanotubes and particularly to discuss 
concrete physical conclusions. This discussion requires a special publication (see |HS1) and 
probably not a single one. Here we only want very briefly to outline the structure of semi- 
classical asymptotics for some natural problems. 

Wave trains propagation. It is natural to formulate the problem of the wave train 
propagation for the reduced equation (15. b|) . We consider the general case when £g and are 
matrices. Let be the certain eigenvalue of the matrix Cq. If is degenerate, to restrict 
on the subspace corresponding to it, we can apply the reduction of §3 (by parameter h instead 
of /i) once again. So we have to restrict the operator to the subspace corresponding to the 
eigenvalue After the reduction our problem has the following form: the leading part of 
the matrix Hamiltonian is proportional to unity r' x r'-matrix with coefficient Hj^g, where r' 
is the multiplicity of degeneracy of eigenvalue , and is a r' x r'-matrix. Thus we always 
can reduce the problem of the construction of semiclassical solutions to the problem with 
the leading term proportional to unity matrix. We construct asymptotic solution assuming 
that the initial problem is reduced to this form. In order to not overload the notation we 
will write for the matrix restricted to the proper subspace corresponding to eigenvalue 

Ijh 

HeS- 

Let ^jJ^{x,0) = 'ipo{x) = exp ^^^ 7 ^) '^o{x), where the phase So{x) is a smooth function 

and (po{x) is 2 r-dimensional smooth vector function of a; G M with compact support M.. 
Then in the semiclassical approximation the solution to dEHI) is determined by solution to 
the Hamiltonian system (j5.7jl with Hamiltonian where depending on 

regimes of longitudinal motion are dehned by formulas (|EHD,(ISinD with initial conditions 
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p^(0) = dSo/dx^Xo), a;(0) = Xq. Denote its solution = 7r(t, xq), x = ^{t^xo). Suppose that 
J{xo,t) = ^ 7 ^ 0 for t G [0,t*] and xq G A4. Then the equation x = ^(xo,t) is uniquely 
solvable, xq = Xo(a;, t) G Af, and, for t G [0,t*], the asymptotic solution to the Cauchy 
problem for system (EH) with initial condition ' 0 |i=o = h) has the form: 


'ip{x,t) = 


exp 


iS{xo,t)\ I ^p{xo,t) 


h ) \^yj{xo,t) 


+ 0(h) 


S{xo,t) = S'o(xo) + / ( X 




dwi 


xo=Xo{x,t) 

HesiT^,^,t)] dt. 


Function ip{xo,t) satisfies the following equation: 

d(p 


dt 


+ iCi^p = 0 


(5.13) 

(5.14) 


(5.15) 


with initial condition (p(xo,0) = <po(3^o)- 

At t > t* the asymptotics of the solution is specihed with the use of the Maslov canonical 
operator (see IHBlEnj ) on the curves Aj = {x = ^{xo,t),p = 7r(xo,t)}: ip = Ka^'^o- 

Remark. In the case of a finite effective potential the effect of splitting of incoming 
wave train into two space components can appear. They are partially reflected and partially 
transmitted, containing harmonics E < max and and E > max respectively. 

Plane wave scattering. Consider the nanotube having the following structure. For 
X < x_ and x > x^ (x± = const), this tube has rectilinear axis, constant torsion angles 
and <F+, and the contraction coefficients D_ and D+ at the ends. Suppose that H = 0, Uext 
and Uint do not depend on t, next(R'(3^)) = ^ ^ > ^+}) Ant{x,y) = 

{Dnt(2/) ^ Dnt(l/) ^ > ^+}- Then for x < x_ and E > + e'P, the system has 

a solution of the form exp{{—iEt + ip_x)/h)[y), k = 1,..., 2r, which represents the plane 
wave propagating along the tube axis with vector (spinor) amplitude g^^{y)- Vector g'P^{y) 
and e'P are the eigenfunction and eigenvalue of the following problem, correspondingly: 


g’P^iy) = e'Pg'PPp{y), no = (5.16) 

M = aa,R(^(9inint)52 - {d2Vmt)d^ + ni(a 2 nint)p - n2(5inint)p. (5.17) 


Remark. Note that eigenfunction g^^ is not the product Wj{y)®p of “pure state” inside 
the term Wj{y) and “pure spin state” pk (hfc doesn’t depend on y). It means that we can’t 
separate spin states and states inside one term. Since the parameter a -C 1, to construct 
guk ggg perturbation theory. To construct solutions in the case when Uint = yl + y^ 

(parabolic confinement) one has to use second order of perturbation theory. 

As exact solutions “at infinity” are not products of “pure states,” we need to expand 
them with respect to the basis Wj{y) ® . We obtain: 

g'P^{y) = \\w\{y),...,wp{y)\\®Es^p~, (p“ = {(p7,..., (5.18) 


The evolution of “initial scattering data” for reduced equation is determined by the following 
transport equation: 


(f ^ (f , t —»• —oo. 


^ -ry _ n 

— + iC,ip - 0 , 
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(5.19) 









The final scattering data for reduced equation is = linii^oo So we obtain the 

final scattering data for original problem in the form: 

gf(v) = \\wX(y),...,w''Ay)\\»E.^*- ( 5 . 20 ) 


In the semiclassical approximation this problem has 2r-dimensional family of asymptotic 
solutions of the form E) ~ Kai(£;){‘^(^)}, where Kai(£;) is the Maslov canonical operator 
on the nonclosed curve A^(£') = = E : x = a;(f)}, t is the proper time 

(parameter on A^(ii^)): dx/dt = ±^y2{E — v^ff{x)). Let E > v!^ff{x) for each x. Then with 
an accuracy exponential with respect to h —0, we have a passage of the incident wave above 
the barrier; as a: —> ±cx) 


'ipix, E, h) -^= exp (j^P±xj P± = \J2{E e''_). (5.21) 

If < max (a;) at some points of the tube axis, then the incident wave is reflected 
off the barrier with an accuracy exponential with respect to h. In the domain x < Xf{E), 
where xj is the rotation point at the energy level E = v^g^Xf^E)), tfj^x^E^h) is the sum 
tfj- + E,h) of the incident and reflected waves 'ip±{x^E,h) with x < Xf{E) (at 

X > Xf{E) 'ijj{x, E, h) = 0{h°°)). As X —>■ —oo, we have 




Vp 


exp 


±^pa; 




± 


(5.22) 


The form of the nanotube after barrier modulo exponentially small corrections doesn’t in¬ 
fluence the solution. Therefore the part of the tube behind the barrier can be removed. One 
can see from the formulas (lO^ that the barrier appears not only because of the external 
potential but also because of the narrowing of the tube {D{x) 0). Thus the constructed 

asymptotics model the situation when the end of the tube narrows conically-like, this result 
in the appearance of a barrier, that is a turning point in the system (15. bp . It results in sharp 
increasing of the wavefunction in the neighborhood of the conical end of the tube, this is 
probably related to the effect of the luminosity of the tube end. 

Bound states. Asymptotic eigenvalues. Suppose that the potential Uext and the 
magnetic field H do not depend on t and the effective potential ngg(R(x)) has a stable 
minimum point Xq. In its neighborhood ngg(R(a;)) has the form of a potential well, which 
generates a family of T(ii^)-periodic trajectories E), 7i(t, E) of system (j5.7|l parametrized 
by the energy E = -|- (a:). Substituting them into we obtain the system of the 

form ip = —G{t)p> with unitary matrix T-periodic with respect to t. We can form a basis 
in the space of their solutions from vector-functions of the form z^{t, E) exp{if5^{E)t), where 
j = 1, 2,... 2r, z{t, E) are T-periodic in t, and the Floquet exponents P^{E) are real for all 
E. We choose them in such a way that \f5^ are minimal. Let be the levels determined 
by the Bohr-Sommerfeld quantization condition 


^min 

^^ 7 id^=^j ^ 2{E - v^^)dx = h + 0 , n = 0,1,2,... (5.23) 

^min 

Then the numbers TJ” = -|- hj3^{E'^'^), where j = 1,2,..., 2r, give spectral series 

(sets of bound states) of the operator £; the 0(h^)-neighborhood of TJ”(h) necessarily 
contains a point of its spectrum (continuous or discrete). Namely, if the spectrum of the 
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original problem on the interval Eq — e^Eq — e) is discrete, then the nnmbers E'^'^ give the 
asymptotics of some of its eigenvalnes. If the spectrnm on this interval is continnous then 
these nnmbers apparently approximate the exponentially small bands of the spectrum (c.f. 

EH). 

The eigenfunctions corresponding to Ej^{h) are determined with the use of the Maslov 
canonical operator. 

5.4 Remark on rigorous justification of constructed asymptotic 
solutions. 

It is natural to discuss the important question about the strict justification of asymptotic 
solutions, which can be constructed by using formal procedures suggested above. In this pa¬ 
per, we almost do not touch upon this problem (see §5.1) whose solution in general situation 
is not trivial and requires an additional investigation. To study this problem one can use at 
least two possible ways. 

1 ) One has to prove that the asymptotic solutions are close to the exact solution of the 
initial problem under some conditions on the coefficients of the initial equation and on the 
function classes to which the solutions of the reduced equation belong. This way is based 
on the technique of obtaining a series of estimates with respect to the parameters fi and h 
from different diapasons. 

2) The second way of justification is to estimate the accuracy of functions obtained, which 
approximate the exact solution of the initial equation. This way seems to be more useful, at 
least from the pragmatic viewpoint, since the obtained explicit asymptotic formulas for the 
solutions of real physical problems can be used. 

5.5 Several effects in nanotubes generated by their geometry and 
external fields 

Finally let us described very briefly several effects which one can obtain from the analysis of 
the asymptotic solutions constructed above. (Some of them have been already discussed in 

§5.1). 

The possibility to model effective one-dimensional potentials. First let us men¬ 
tion elementary, but curious, properties of the above equations, which are caused by the 
possibilities of the nanotechnology: changing the geometry of a tube placed in a homoge¬ 
neous electric field, one can model different one-dimensional effective potentials. 

First, we consider a tube of constant section whose axis is a plane curve on the plane 
(ri,r 2 ). Let the electric held of strength F^ext be directed along the 0x2 axis. Then the 
effective potential has the form (p = next(R-()5:^ )) = Ee,ctR2{x). If the tube is curved not too 
much with respect to the axis ri, then we have a: ~ ri. Thus, choosing the tube axis as an 
appropriate curve r 2 = r 2 (ri), we can model the potential, the double potential well, etc. 

As an example of a nonplanar tube, we consider the spiral 
^{x) = (Pi = const, p 2 = const 

are parameters) in the held F^ext(0, sina, cosa). The ehective potential contains the 
oscillating and linearly increasing components: ip^x) = — sin aF^extPi sin(a;/'ypf+^) — 
X cos aEextP2/ \/ Pi + P2- If « = 7r/2, i.e., the held is directed perpendicular to the tube 
axis, then we obtain a periodic potential and the equations obtained above coincide in the 
hrst approximation with the Mathieu equation. If a = 0, i.e., the held is directed along the 
tube, we obtain the Airy equation. A more complicated example is the case in which the 
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tube axis is a winding of a torus: in this case, in particular, we can obtain almost periodic 
potentials. Similar results can be obtained by changing the thickness of the tube along the 
tube axis. 

Dependence of the effective one-dimensional potential on the tube thickness. 

It is easy to show that for the above potentials modelling the “soft” and “rigid” walls, the 
effective one-dimensional potential of the longitudinal motion depends on the tube thickness 
do{x) as l/do{x)^. This dependence is a purely quantum effect and is caused by the effect of 
the size quantization. Thus, narrowing the tube, we obtain a barrier in the one-dimensional 
motion, while expanding the tube, we obtain a potential well or a “trap.” 

Semiclassical splitting of the adiabatic term. The “adiabatic” terms can split in 
the semiclassical approximation due to spin and the magnetic held. In particular, if the 
splitting is caused by spin effects, then spin affects the “classical” trajectory. For the case in 
which the adiabatic term is nondegenerate, we have A = 0. If degeneration exists, then the 
momentum matrix A is nonzero. This results in the appearance of an effective “dipole” that 
interacts with the projection of the magnetic held on the tube axis, and thus an additional 
phase of the wave function (the Berry phase) appears. 

An increase in the electron density near the endpoint of the nanotube caused 
by reflection. Since the longitudinal energy decreases with decreasing cross-section of the 
nanotube, in a nanotube with a “closed” endpoint, the wave packet reflects from the closed 
endpoint. In this case, near the endpoint of the nanotube, there is a sharp increase in the 
electron density, which, apparently, is related to the ehect of luminous emittance observed 
in nanotubes. 

The Berry phase of the wave function. The term [<F 3 ;A — (^2111 — Yin2,H)] 0 
EsP^ in the Hamiltonian (EUl) proportional to the momentum can be excluded from this 
Hamiltonian by a change of the wave function. This change has the form 




exp 


Xj{x)dx ) 


(5.24) 


where Xj{x) is an eigenvalue of the matrix $ 3 ,A — (T 2 ni — Yin 2 , H). In fact, the existence of 
this term results in the appearance of an additional phase of the wave function. This phase 
must be consistent with the boundary conditions at the endpoints of the tube, which gives 
a correction to the quantization condition. For example, we consider the Born-Karman 
boundary conditions that are equivalent to the case of a closed tube (the nanoannulus). 
While constructing the eigenfunctions in this case, it is necessary to require that the total 
phase of the wave function be 27r-periodic, with the above adiabatic correction taken into 
account. It is well known that the adiabatic phase responsible for the correction to the 
quantization condition is the so-called Berry phase. 

If the magnetic held is zero, then the Berry phase is reduced to exp Aj{x)d^'j, where 
Aj{x) is an eigenvalue of the momentum matrix A. 

Ultrashort modes and the state density. The correction to pxi becomes comparable 
with 1 in the case h = This means that the adiabatic asymptotics of the form (1071) 
ceases to hold. However, for h < p and for the same total energy, there exists a state of 
the form (1071) that belongs to the next subregion of the transverse quantization. From 
the physical standpoint, the fact that there is no asymptotic solution of the form (1071) in 
the case of ultrashort waves in a curved tube means that such fast modes are destroyed in 
passing from a straight tube to a curved tube; they “fall” to the next subregions because of 
the interaction with the tube “walls.” 

The “instability” of ultrashort modes with parameter h < in a curved tube can lead 
to a change in the density of states and, in contrast to the straight nanotube, ehectively 
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increase the Fermi level (cf. j92ji. 

Spin beats. The operator Lay in the semiclassical Hamiltonian, which corresponds to 
the spin-orbit interaction, leads to splitting of the adiabatic term into 2r semiclassical terms 
of the longitudinal motion. To each of the semiclassical terms, there corresponds, in general, 
its own phase Sj{x,t), j = 1,... 2r. If the variation in the phase due to the spin-orbit 
interaction is small, then the phase can be expanded in a series in the constant of the spin- 
orbit interaction a. The zero term in the expansion corresponds to the classical motion of 
the “spinless” particle, and the correction results in the appearance of an envelope for the 
fast-oscillating exponential. Apparently, the appearance of such an envelope is related to 
the observed effect of the electron density beating along the tube 1231' 

6 Concluding remarks 

Let us briefly summarize the results of this paper, sometimes repeating the above arguments. 
The adiabatic approximation is one of the main tools for analyzing the solutions of linear 
stationary and nonstationary problems in modern mathematical and theoretical physics. 
Different versions of adiabatic approximation are used in problems of molecular physics, 
solid-state physics, plasma physics, theory of internal and surface waves in fluids, averaging 
theory, quantum gravitation theory, etc. 

The adiabatic approximation is used in situations in which the study of some classes of 
physically interesting wave process described by a “large” system with N degrees of freedom 
can be reduced, on some time intervals, to the study of a simplihed “effective” system with 
n < N degrees of freedom. In this case, on the one hand, a sufficiently wide range of 
states and solutions with some general characteristic properties (but not a set of individual 
concrete states or solutions) is considered, and, on the other hand, it is not assumed that 
all the solutions of the original “large” system can be described in the same way. From 
the physical viewpoint, this possibility is usually ensured by the fact that the problem has 
different spatial or spatio-temporal scales, which, from the mathematical viewpoint, means 
that there is a small parameter characterizing the different scales of the problem. As a rule, 
the different scales are manifested in different dependences of the coefficients of the original 
equation or boundary conditions on various variables (or groups of variables), i.e., on {N — n) 
“fast” and n “slow” variables. Thus, it is natural to divide the study of such distinguished 
wave processes into two stages: 1) to derive “effective” reduced systems, 2) to hud their 
concrete solutions and then to reconstruct the total solution describing the entire process. Of 
course, these considerations appear in different helds of physics and mechanics; the problem 
is to realize them in concrete formulas. Moreover, in each held, its own terminology is used. 
For example, such reduced systems correspond to terms in molecular physics, to modes 
in hydrodynamic problems, to subbands of dimensional quantization in nanophysics, etc. 
For various reasons, it is convenient for us to use the terminology accepted in physics of 
low-dimensional systems. 

In turn, roughly speaking, the idea to derive simplihed systems has two stages: 1) hrst, 
after the n above-mentioned degrees of freedom are “frozen” (it is assumed that the diher- 
entiation operators with respect to the slow variables commute with the slow variables), the 
operator determining the ehective equation can be obtained as an eigenvalue of an auxil¬ 
iary spectral problem with N — n degrees of freedom; realizing this stage, we otbain ehective 
Hamiltonians (or dispersion relations), well-known in molecular physics; 2) the operator (the 
“quantization” or the Peierls substitution) determining the ehective reduced equation with n 
degrees of freedom is reconstructed simply by taking into account the fact that the diher- 
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entiation operator with respect to slow variables does not actually commute with the slow 
variables. We note that the distinguishing of the “frozen” (slow) variables can be natural and 
obvious as, for example, in problems of molecular physics or in mechanical problems with 
different spatial scales, but can be more “latent” as in electron waves in crystals and in the 
averaging theory. In this case, an additional degree of freedom appears in the regularization 
of the problem. 

In some problems, it is sufficient to perform the “naive” quantization of the effective 
Hamiltonian (of the dispersion relation), but, in many problems, a more accurate analysis 
is required and certain difficulties arise (see, e.g., 0, §56). We note that, as a rule, it is 
impossible to write the effective reduced equation exactly; the problem is to construct a 
minimally reasonable number of terms in the expansion of the operator in the effective re¬ 
duced system, which allows a correct construction of the equation describing a wide range 
of wavelengths. Indeed, the most popular approach in the adiabatic approximation, orig¬ 
inating from Born and Oppenheimer works, is based on the assumption that the desired 
solution depends smoothly on all the variables. Thus, for example, only the lower energy 
levels are usually “grasped” in spectral problems. At the same time, in many physical situ¬ 
ations, for example, in describing the valence electrons in crystals, the higher energy levels 
are most interesting, but, strictly speaking, they cannot be considered under this approach. 
In this case, one can use the semiclassical theory proposed by Maslov, which, however, is 
based on the assumption that there is a sufficiently rigid relation between the excitation 
level and the parameter characterizing the different scales of the problem. From the math¬ 
ematical viewpoint, this means that, in the problem under study, along with the parameter 
characterizing the different scales, which is naturally called an adiabatic parameter, there 
is another (“semiclassical”) parameter characterizing the excitation level of the state under 
study; in this case, the form of the approximate (asymptotic) solution depends, as a rule, on 
the relations between these parameters. This can easily be verified comparing the averaging 
method, the Born-Oppenheimer method, and the Maslov method, which give solutions with 
“slow,” “medium,” and “fast” variations (with respect to the slow variables). 

In this paper, we propose an approach based on the above considerations. This approach 
allows one, hrst, to describe all the states listed above, in the range from slowly varying 
states to fast or, sometimes, even “superfast” varying states, and, second, to classify them 
appropriately. In particular, this approach explains why the states obtained by the Born- 
Oppenheimer and Maslov methods can be treated as the states on the same subregions 
(terms) corresponding to different “longitudinal states.” So the effective reduced equation 
thus obtained describes not only the states at the “bottom” of the subregion, but also the 
states corresponding to the higher energy levels. We note that, without this reduction, 
choosing the subregion for higher energy levels can be a rather complicated problem by 
itself. 

The approach proposed here is based on Maslov’s observation that the problems in which 
the adiabatic approximation is used can be interpreted as problems with operator-valued 
symbol, and these problems are well known in mathematical physics. This means that the 
operator determining the original “large” system is a function of two groups of operators with 
“small” and “large” commutators (or “large” anticommutators, as in the electron-phonon 
interaction problem). The use of the concept proposed in this paper, which is based on 
equations with operator-valued symbol, allows one to deal with different linear adiabatic 
problems from the unified viewpoint. 

The realization of our approach is based on representing the solution in the form (3.5) 
and obtaining the effective reduced equation (3.7), which is a generalization of the Peierls 
substitution. These formulas, together with the algorithm described in §3, give one of 
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the main results of the present paper. In these constructions, the key role is played by 
the techniques of noncommuting operators based on elementary notions from the Maslov 
operator method. The constructed algorithm allows us, with any prescribed accuracy, to 
calculate the operator determining the effective reduced equation (3.7) and the intertwining 
operator reconstructing the solution of the original problem from the solution of Eq. (3.7). 
In particular, the algorithm thus constructed allows one to obtain accurate estimates for 
the minimal number of terms necessary to construct the leading part of the asymptotic 
solution. It should be noted that, in the reduction procedure, the possible degeneration of 
the term (the effective Hamiltonian) must be taken into account. In §§4-5, these formulas 
are illustrated by several examples from different helds of physics and mathematics. Some 
heuristic arguments leading to formulas (3.5) and (3.7) are given in §2. 

In §3, we also show that the “operator separation” of variables can be treated as a 
“quantum” (or wave) analog of the procedure of excluding (holonomic) constraints in classical 
mechanics (isaisii and others). Indeed, the imposed quantum constraints can be treated 
as the restrictions arising due to the conhnement potential in the ambient conhguration 
space. In this case, it is natural to assume that the “condition of dimensional quantization” 
is satished, i.e., the wavelength in the directions normal to the manifold corresponding to 
the degrees of freedom of the effective system (i.e., the “limit” manifold) is compatible with 
the “width of the him” surrounding the limit manifold. Here the most important is the fact 
justihed in this paper that, excluding the constraints, one can, in general, obtain different 
effective Hamiltonians depending on the energy of the “longitudinal” motion. In the case of 
fast oscillating longitudinal states, this leads to different classical Hamiltonians determining 
the motion on the “limit manifold.” 

The study of solutions of reduced equations on some distinguished subregions is the sec¬ 
ond part of our approach. Depending on the relations between the parameters, this study 
can generally be performed by different methods. In §5, taking into account the fact that the 
problem contains two parameters, an adiabatic and a semiclassical, we classify the solutions, 
and this classihcation shows that the excited states are constructed differently than the lower 
states. In the construction of excited states, the momentum in the intertwining operator 
cannot be neglected, i.e., the intertwining operator is indeed an operator and cannot be 
replaced by a function, as it is usually done in different versions of the Born-Oppenheimer 
method (in particular, in solid-state physics). The existence of the momentum operator in 
the intertwining operator shows that, from the viewpoint of the Born-Oppenheimer method 
for excited states, a “distortion” of the term occurs. For fast varying solutions, the main 
methods in this case are the semiclassical approximation and the WKB-method; if there are 
turning points and caustics, then the WKB-Maslov method is used. It is well known that, 
in the realization of this method, one must pass to classical Hamiltonian systems. One of 
the elementary, but important, consequences is the fact that the classical systems can be 
different for different excitation levels and the “small” terms in the original equation can sig- 
nihcantly affect the semiclassical characteristics for some values of the longitudinal energies. 
In particular, in the case of a degenerate adiabatic effective Hamiltonian, the degeneration 
can be removed in the semiclassical approximation; in this case, the adiabatic term “splits” 
into several semiclassical terms (Hamiltonians). We consider an example (nanotubes) and 
show that the interaction of spin with the conhnement potential can change the classical 
trajectories of the longitudinal motion. 

As was already noted, the possibility to obtain numerical solutions, graphs, etc. at this 
stage signihcantly depends on the specihc character of a concrete problem and requires sep¬ 
arate publications. In the present paper, we briehy describe this procedure for problems 
related to the modern held of nanophysics and restrict ourselves to rather general formulas. 
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In §5, we derive some simplest conclusions for models arising in nanophysics. In particular, 
we show that, placing nanotubes of various geometry in a constant electric held, we can 
model various one-dimensional potentials, for example, “double well” type potentials, peri¬ 
odic potentials, etc.; the degenerate adiabatic term (for example, in the case of a tube of 
circular section) can split into several semiclassical terms (effective Hamiltonians), etc. The 
problems concerning detailed studies of how the spin affects the classical trajectories, the 
electron density pulsation due to spin, etc. are beyond the framework of this paper. 

We believe that the arguments and formulas given in this paper can be helpful in studying 
the problems arising in solid-state physics, hydrodynamics (waves on water), the theory of 
shells, plates, and rods, and in nanophysics. It seems possible that this method can be used 
in weakly nonlinear situations. 

In conclusion, we make a remark concerning the list of references. As was already noted, 
the number of works dealing with the adiabatic approximation and its applications is ex¬ 
tremely large; our list of references does not absolutely pretend to be complete; here we 
present only several papers that are to some extent close to our approach. 
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